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(ii) The species scattering from a centre is supposed to distribute itself 
uniformly in all directions. The average distance through which an individual of 
the species moves from habitat to habitat will be spoken of as a “flight,” and there 
may be such “flights” from locus of origin to breeding ground, or again from 
breeding ground to breeding ‘ground, if the species reproduces more than once. 
A flight is to be distinguished from a “‘flitter,” a mere two and fro motion associated 
with the quest for food or mate in the neighbourhood of the habitat. 


(iii) Now taking a centre, reduced in the idealised system to a point, what 
would be the distribution after n random flights of N individuals departing from 
this centre? This is the first problem. I will call it the Fundamental Problem of 
Random Migration. 


(iv) Supposing the first problem solved, we have now to distribute such points 
over an area bounded by any contour, and mark the distribution on both sides 
of the contour after any number of breeding seasons. The shape of the contour and 
the number of seasons dealt with provide a series of problems which may be spoken 
of as Secondary Problems of Migration. 

A little consideration of the Fundamental Problem showed me that it presented 
considerable analytical difficulties, and I was by no means clear that the series of 
hypotheses adopted would be sufficiently close to the natural conditions of any 
species to repay the labour involved in the investigation. At this stage the matter 
rested, until last year Major Ross put before me the same problem as being of 
essential importance for the infiltration of mosquitoes into cleared areas, and asked 
me if I could not provide the statistical solution of it. He considered that we 
might treat a district as approximately “equi-swampous,” and thus my conditions 
(i), (ii) above could be applied to obtain at any rate a first approximation to 
the solution. 

Starting on the problem again I obtained the solution for the distribution after 
two flights, an integral expressing the distribution after three flights, which I 
carelessly failed to see could be at once reduced to an elliptic integral, and the 
general functional relation between the distribution after successive flights. At this 
point I failed to make further progress, and under the heading of ‘‘The Problem 
of the Random Walk” asked for the aid of fellow-mathematicians in Nature*. The 
reply to my appeal was threefold. Mr Geoffrey T. Bennett sent me in terms of 
elliptic integrals the solution for three flights. Lord Rayleigh drew my attention 
to the fact that the “problem of the random walk” where the number of flights 
is very great becomes identical with a problem in the combination of sound ampli- 
tudes in the case of notes of the same period, which he has dealt with in several 
paperst. Lastly Professor J. C. Kluyver presented a paper to the Royal Academy 

* July 27th, 1905. 
+ Phil. Mag., August, 1880, p. 75; February, 1899, p. 246. 
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of Sciences of Amsterdam, entitled ‘‘A local probability problem.”* Professor 
Kluyver obtains the general solution in terms of the integral of a product of 
Bessel’s functions of the zero and first orders. He deduces Rayleigh’s solution for 
n large, he shows that the Bessel function integral represents a series of different 
analytic functions in different intervals, and proves a number of special problems 
of very considerable interest. Referring to his general solution, he writes, however : 


“From this result we infer that the probability sought for is of a rather intricate 
character. The n+1 functions J are oscillating functions, and have their signs 
altering in an irregular manner as the variable wu increases. Hence even an 
approximation of the integral is not easily found, and as a solution of Pearson’s 
problem it is little apt to meet the requirements of the proposer.” t 


Kluyver’s solution is of extreme suggestiveness for the analytical theory of 
discontinuous functions. In the endeavour to express it in a form suited to my 
special purposes I have come across a long series of Bessel function properties, 
some at least of which seem to me novel, but have unfortunately no bearing 
on the problem of migration. If we turn to Rayleigh’s solution for n large, I 
must confess at once to being unconvinced of the adequacy of the proofs used to 
deduce it, especially that in the Theory of Sound. Kluyver’s proof of Rayleigh’s 
solution § appears to me to also require much strengthening, and in neither case do 
we have any practical measure of what the number of flights must be before we 
have in practice a reasonable accordance between the discontinuous Bessel’s function 
integral expression and the Rayleigh solution of Gaussian frequency type. 

After a good many failures I have succeeded in obtaining a solution in series 
of the Bessel function integral, but this not of a character to be of service for 
Srequent arithmetical calculations. It serves, however, to test the approximation 
of the Rayleigh solution and the accuracy of the solutions for few flights obtained 
by other processes. At this stage I realised that the functional equation between 
the distributions for n and 7 + 1 flights could be solved graphically, and that starting 
with the known distributions for n= 2 or n= 38, we could by very great labour, 
but absolutely straightforward graphical work and the use of mechanical integrators, 
build up in succession the solutions for n = 4, 5, 6, 7, etc. I proposed that this 
process should be continued until the graphically found distribution coincided with 
the plotted values obtained from the above solution in series. This was achieved 
for n=7. For n=6 and n=7, the solution in series approaches to the Rayleigh 
solution, with which for all practical purposes it may be asserted to cvincide for 
n=10. We have thus reached a continuous graphical illustration of the transition 
of a series of discontinuous and, in many respects, remarkable analytical functions, 
step by step with the increase of 2 into a normal curve of errors. The relation- 


* Koninklijke Akademie van Vetenschappen te Amsterdam. Proceedings, Oct. 25, 1905, pp. 341—50. 
Tt loc. cit. p. 348. { 2nd Edition, § 42a. § Kluyver, /oc. cit. p. 345. 
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ship is a noteworthy one, and not without suggestion for other branches of 
investigation. 

The exact method of graphical solution will be described later ; the whole labour of 
it, involving many weeks’ work, was due to my assistant, Mr John Blakeman, M.Sc. 


(2) . General Analytical Solution of the Fundamental Problem. Let the origin 
be taken at the centre of dispersion and r be the distance of any small elementary 
area a from the centre of dispersion. Let ¢,,(7°). a be the frequency of individuals on 
a after the nth flight, and ¢,,,(7°)a their frequency on the same element after 
the (n+1)th flight. Let / be the length of the flight. Then only those individuals 
who were on a circle of radius / round the centre of a after the nth flight can reach a 
with the (n+1)th flight, and only those individuals of these who take their flight 
in one definite direction. Let O be the centre of dispersion, C’ the centre of a, 
P a point on the circle of radius / round C, 4 PCO=86, then the frequency per 
unit area at P is ¢,(7°+/—2rlcos@), and the amount which goes in directions 
between # and +60 is d@/27. Hence the frequency per unit area at C’ after the 
(n+1)th flight is given by: 


Ae a | "dh, (ASP EIOT] con O\ des mee ae eee (i). 
277} 9 ; 


This is the equation, which I shall speak of as the general functional relation between 

the densities at successive flights. Now assume: ¢, (7°) =C, J, (uv), where C;, is any 

undetermined function of n, J and wu, and w is at present an undetermined variable. 
Then by Neumann’s Theorem *: 


J, (uv? +P — 2rl cos 6) =J,(ur) J, (ul) +2 S J, (ur) J, (ul) cos t6. 
1 


Hence : x |, J, (uJ +P —2rl cos 6) d0=C,, J, (ur) J, (ul) 

Sas Ig 1) 2 Dye 
Therefore C= J, (ul) Ce 
It follows that C,=D {J, (wl)}", where D may be any function of /, but not of n. 
Thus we have: P,P) = DiI (Und (uly 


where we may sum for any or all values of w. 

Now when n= 1, ¢, (7°) must be zero, for all values of 7 except r=/ to /+7, and it 
then equals N/(2z/7), + being very small and MN the total number issuing from 
the centre of dispersion. We know, however, thatt: 


| au] upf (p) Jn (up) Jn (ur) dp=f(r), if g<r<p; 
0 
4 = 0)" \ifae> p08 <ig. 


* ©. Neumann, Theorie der Besselschen Functionen, 8. 65. 
+ Gray and Mathews, 7'reatise on Bessel’s Functions, p. 80. 
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N 
Now take n=0, g=l, p=l+r and toa ek 
then we have : une J, (ul) J, (ur) rdu = @, (7°) 
; 0 Qelr i ‘ ; 
or o (pe a l) J, (ur) d ii 
A aos pe NOT EVIIT S| hein aire erie Seam Sod): 


This determines the form of D and the summation of wu; for, if we take 


by(o?) =e | ase (Csa Ne CUNC)) eenemecteeercreoes ) 


we satisfy the general functional relation (i) and further the initial equation (ii). 
Let P,, (7) be the probability that an individual after n flights will be a distance r 
. or less from the centre of dispersion. Then clearly 


P, (vr) =2 | td baa (1) 
0 


=¥|'r ar | uJ, (ur) {J, (ul)}” du. 
0 0 


st es Nea 
Hence Game au |” um yeu (ur) ur} {J, (ul) 
d (ur) U 


a ni. rJ, (ur) {J, (ul) "du, 


Orv =U: =v] J, (v (2) Cielo un a (iv). 


(iv) is Kluyver’s fundamental solution, which he reaches by a very different 
and more general analysis. (iii) is the form of it which best suits my present 
investigation. 


(8) On an expansion in series of the expression for >, (7°). By straight- 
forward but somewhat laborious multiplication it can be shown that: 
= 3, (6n—11 
{Jo (2Vy) et} = 1 — pny? — pny’ + oe 
(50n—57)m,,_ (1892 ~21250+270n") n 


T8004, op 103,680 ae 


* Gray and Mathews, loc. cit. p. 13. 
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Hence putting 2 Jy =2, 
‘fh een ee Ol ee 


n__ ,— nz mea ag A 
le: Re EG 192 256 
(50n—57)n 2° — (1892—2125n+270n*)n 2” . 
1800 1024 103,680 4096 7 tf (v) 


Sees {1 -—az—a,2° — a2 — 4,2" — 4,2" — ete}, 
let us write, for brevity. The a’s are then known coefficients. 
Now by (iii) 
1 D 
one | uJ, (ur) {J, (ul) }" du. 


0 


But we know that*: 


@ ; ‘ ‘ i 

iL ue? 7 (ur) du = -. Pie eae RRR, one ne (vi) 

Write: DE eat amare nC rears Cut aoe ae eens (vil). 
@ 

Thus : | ue 2”? J, (ur) du= * io” AT ney ee ee (viii) 
0 


Differentiate (vill) s times with regard to o°: 


; 2s+1 ,—$u?0? ds Ih eye 
(-4) Cr Jy(ur) du= Frey {ase 2 a 


0 


Hence, if B= —20°/r", we have: 
@ 
thee 
1=| wath e 2” " 7 (ur) du 
0 


sa eee 
> dB B 
O21 : 

~ eats Vos) say. 


We have therefore: 
N ert UNIS Che OST 


Pn (7°) “a ono? Qa res Ms 


ge ; Sear #) 
where it remains to evaluate tes = 9 Bi ( a ; 


We find: 


: 1 /r\* —72/9¢2 Fp th WN 
4,= —=(—Je 2-2—4+-— 
: 8 (“) : o” 8 Av 
1 /P\*. 72/902 Oa ee “) 
= ne 6-9 — Sie 
) | : o = 46° So 
4 7 l 7 


a\ 10 ss 72/202 ( 7? r 
2 DA de IS eee ae 
ye ‘ Be: ies Suis 


i 
| 
| 


S 

i 

| 
| 
a= 


* Gray and Mathews, loc. cit. Eqn. (162), p, 78 


A MATHEMATICAL THEORY OF RANDOM MIGRATION 9 


5 1 Pann \ Aa — 72/26" { a2 at 6 9 8 | 4l0 
Ane () e *Pe* (120 — 800 5 + 150 “i -25 en Cat a ), 
oO 


«64 \o 160° 320° 
1 /r\* -19/202 / ie re ei (eI fr \? 
~= —-- (— Oe 16019505 — +2 : 
te 128 (=) ; (7 v ae fee o* ce o° 4 8 -o4" 8 (=) = 64 (=) ) 


Remembering that by (vii) /?=20°/n, we have from (ix) 


. Ne kala? 1 raw lene 1 ip OSES ae ed 
ee eae [Teen me = = = 
$n (7) ono An ( z o a? 4 ) on? (6 ? o a 4a° 8 =) 


6n—11 ( r a 7 1 7 
eh 4 Vogt) ee 
Gi 1 Aa : en : o ao zu 16 =) 
50n — 57 he an (i PAS ee 
———_(]20 — —+150——-2 : Sant 
a 1800n+* ( ae ee mM o* 3 oa x 160° 32 =) 
1892 —2125n+ =e 7” 4 r 229578 9 pi 1 = 
iooeitee ee Gig GF B® G4g™ 


\ 


This is the general expansion for the distribution of the individuals emerging 
from a centre of dispersion after n random flights. Clearly if we want to go as 


1 ° ¢ 2\9 ° 
far as — we must retain terms up to (7°/o*)", and the convergence is small for n 
n 


small. Thus for the first two or three flights, (x) as far as I have calculated the 
terms gives poor results, even if they are notwithstanding better than the Rayleigh 
solution. The arithmetical work required to calculate the ordinates is also severe. 
If we put n=, we have 


Lord Rayleigh’s expression. Now o*=4n/’, hence unless / becomes indefinitely small 
as 2 becomes indefinitely large the population becomes widely scattered. If the 
single flight 7 be very small, but the total flight nl be finite, then }n/* tends to 
become vanishingly small, or the population remains close to the centre of dispersion. 
This is really the ‘‘flitter” as distinct from the flight. 

Examining the solution found it is clear that it may be looked upon as the 
sum of products of two factors, one series of factors not involving r/o but only n 
and the other not involving n but only r/o. Thus we may write 

dn (1) = NV (v1, + ,.0, + V,0, +50; + -.+); 
where 


ili —172/g2 
= € 2 / ; y=l, 


Lp edie? Li ies ; 
ea cae (1-35), v,=0, 


bo 
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b Pea 92 geal re 1 
o,= BAG, 2 |e 2-2 se hee 3 i — 
270" oe AG An 
Ty Reeetiee ei Ae cg wel 1 
On= -e2 oN Ga ,+ j-5 a) y= =, 
Qiao” OG 240g | ee.0, on? 
OG. 155.25 ees hewangiint a. toe PIER ane eee ete pecan a eee te (xii). 


The w-functions form a series of such special interest that a few of their 
remarkable properties will be stated in the next section. 
(4) Properties of the w-functions. 


Let us consider the pth moment round the origin of the 2sth w-function taken 
over all plane space. We will denote it by m,.,. Then 


Qar oO 
nae | dé [, rdra,, 7” 


= 25. | ” cor? dr s Avis Sas Dy eecuecan ge ca ha tincte Cr se (xiii) 
0 

Be ae py ne) 

Now Ga Soa B) dp (ae Reete eee (xiv), 
and by B= —20°/7’ we have dp =" dr. 

Hence writing p=2q we find 

oad 2 : s-qr1 dé /1 

Mg «p= (— 1)? (Dot? i ner ( 2 er) A habe aiid, WiPe ay 


Integrate by parts and we have 


, , eee ds7 1 i Eihe ko) Ps, ds 
= (—1)1*8-(2g? erg=a (pie PAY omer e—q-2 
tga =(— 12009 | {Br Te (Bel) “—(e—g—1)[ Bree 
The part in curled brackets vanishes at the limits and thus 
May (1) (20°)"(8—q=1) [Br ap (ge) 48 


= Meg, 99-9 (S—Y — 1). 


Repeating this process we find 
Mq,2= (8— 1 —g) (S—2—g) (S— 3g)... (— 4) 
x (1) x (20°)! x i A Ur Ben 
-0 


The integral is finite and known; hence if g be less than s we find for integer 
values 


ie) 


Man ge, GA S> sissies ca See eee (xvii). 


Now consider w,, as made up of two parts, 
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Then it is clear that y.,, if g be less than s, consists of powers of 7 less than 
s, and therefore 


| ONO Tar = 0. 
0 


Accordingly a remarkable property holds for the y-function part of the w-function, 
namely, if y,, and x. be two such functions, then it follows that 


| : ey. ya rdr=0, Wogeandag (be diferent, <.5~- 00.5... (xix). 
Returning to (xvi), let us put g=s, then 


Maa= — 20") |s |" Be ehdB 
hes 


2a°)§ |s 2 
~ i ®(-1) I. ao dar, 
or, Hie a 1 a2 ((|15) es. Se sated ec yen ny cence (x 


Let us now consider the integral over the plane 


@ 
L=%F | WagXo5 TAY. 


0 


Except for the last term in y,,, it will consist of a number of terms having for 
factors m,,,, with g<s and these all vanish. The last term in y,, is 


enaee 


Tee on | pages Toi ee | oe 
0 


and accordingly 


Ye o> 


or by (xx) ii 


is) RE MPN a cu oe fe an ine (xxi). 


Hence we have the following properties : 


(a) The integral all over the plane of distribution of one product of a 
x-function into an w-function of a different order is zero. 


(b) The integral all over the plane of distribution of the product of a 
x-function into an w-function of the same order is, if 2s be the order, equal 
to (\s)’ 


These properties enable us—as in the case of Bessel’s functions or Legendre’s 
functions—to expand any function symmetrical round a centre and a function 
only of the square of the distance from that centre in -functions. 


s=0 


Thus let VMs (D355); 


s=0 
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multiply by x,, and integrate all over the plane, 
Qa | y EF (r) x.7dr =b.2ar : O3XnTOr=0;, {\s}’. 
0 
Hence aia Ei? iA setae aaeae et oat XXIl), 
=p}, PO 


Now x5 consists of an algebraic series in (2). Thus the discovery of the 
OF 


value of the integral | Gr (7°) x,,.rdr depends solely on the determination of the 
0 


odd moments of /F'(r*) between 0 and «©. We conclude therefore that an 
expansion in @-functions involves merely the determination of moments, such as 
every statistician has been accustomed for years to calculate. This is not the 
proper occasion to deal fully with the properties of the o-functions, nor to 
generalise them for odd powers of 7, and to consider the convergency of 
expansions in terms of -functions. They will be discussed on another occasion, 
but the present writer believes that they will be found of not inconsiderable 
service, not only in statistical problems, but in certain physical problems where 
intensity round an axis diminishes with the distance. 


(5) Two further problems are of service for the theory of dispersion. 
Suppose 
F(r°)=S Quen) 


s=0 


Integrate over the plane and remember that x,=1, 
2a | “P(r rdr= ‘S 2a | Dos @osX VA 
0 


ben eee ee eee (xxiii). 


Thus the first coefficient is merely the total volume of the surface z=/'(7’), 
taken over the plane. 
Next consider the second moment 


@ s=O (o.0) R 
Qar | PEP) fdr= 8 27 | Car Os. 
0 s=0 0 


Every term of the summation vanishes except for s=0 and s=1, and the left- 
hand side is the second moment of the function about the axis perpendicular 
to the plane through the centre = volume x (swing-radius)’ =}, x A’, say. Thus: 


co) aK 1 2) Syne a2 ‘ 
b,x R= a be BO 8p + al be ate (1 a =| rdr 
or o Jo HA L078 


= 2b,0° + b, (2-4) 0 =2 (b,—D,) o°, 
or y= b, {1 EA eta gam Coenen fe essere case smile stsuin: (xxiv). 
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Thus far no choice has been made of o% If we take o°=1K’, we have b,=0, 
or if o° be taken half the square of the -swing-radius about the axis of the 
solid of revolution z=/'(r’), that is if o be the swing-radius of the solid about 
any plane through its axis, then the second term in the expansion of F (7°) 
in -functions disappears. 

We are now able, I think, to grasp the relation of the Rayleigh solution 
to the complete solution of the random scatter round a centre of dispersion. 
If $, (7°) be the function giving the distribution after n flights, then 4, (7°) 
can be expanded in a series of w-functions, 1.e. 


hy (7°) = b,@) + byw, + Da, + 0. + Dy Gag $ oo 


By choosing the o° of the w-functions=34°, this becomes, since 0b, the 
volume = JN, 


No -0Hf0 
OM eee, 2 eb, +. + Osx, + sf 


Qo 


Lord Rayleigh’s solution provides the first term of this series, or is the 
correct solution to two terms in the expansion by o-functions. It possesses 
the properties (a) that its volume is the same as that of the complete solution, 
and (b) the mean square deviation from the centre of dispersion is the same, 
.€. 20°, as for the complete solution. 

The latter depends upon the fundamental property of the o-functions that 


| Or Or= 0. if s be > I. 
0 


The expansion in w-functions shows us at once that, whatever be the magnitude 
of m, the mean square deviation from the centre of dispersion is Jnl, and this 
gives us readily a rough measure of the range of habitat of any species for 
which n and / are approximately known. 

Another poimt may be noted here as to the Rayleigh solution. That solution 
is the best fitting Gaussian error surface to the distribution, ze. its volume and 
its standard deviation are the same as those of the actual distribution, whatever 
nm may be. If we take the section, however, of the distribution through its 
axis the standard deviation of this according to the Rayleigh solution is o=J dnl, 
but this is not the standard deviation of the section of the actual distribution, 
ae. the Rayleigh solution does not give the best fitting normal curve to the 
section. It gives only the standard deviation corresponding to @,. It is of 
some value to note what are the standard deviations of other component 
w,, terms. 

To obtain this we must determine the area and even moments corresponding 
to any w,, term. Let 


S+5 -3 1 
vee “ond = Sy Ve B ae (ae o) a, 
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whence integrating by parts: 
- +1 By Ve ene | ee ane 
Ay = su(—1) (s 5) (s 5) (8 ya _, 6 ePdB 


im, it IN ff 3) 5 i fp — $9 
= 522 (8-5) (s 5)(s aya € da 


0 
il ={ ;) eho Gi eee 1 
= -(s—— ees aes Se oe AN eR etamtee ee eee (xxv), 
Jono2\ 2 2 Dinan? 
s= or >l 
lis=0 Ag ga (xxvi) 
: (ily hen he ee ea 
I now take: ea | w,,7? dr 
0 


and find, reducing in the same manner, 


1. oe i 3 / "| es 
p= == s—p-—-— —p-—-—)...( —-p+t- he Bis DY ao (HOS I al earl), 
ian peter ( P 5) ( f 7, \ Pro)” Ces ) 


Clearly : 


by (xiii), hence : 


bf as 1 rus 1 
Or, We 4a NOT ‘(s—p-5) (s—p-5)...(—p+5) 


INS 52 Lee ee eee eee (xxviii). 
Thus the odd moments of the o., functions are known“. 


For the particular case when p=1: 


Pips (3-5) (s-2) - ;) (xxxil) 
Pisa. Se 5 Fab Piet fe rare 3 3 


if k,, be the swing-radius round the axis of the function o,,. 


Hence by (xxv) 


a (—4 o aie 
ed = aoe = sy cre (XxXxIii). 


* If x=r/o the following finite difference and differential equations are fundamental in the 
theory of the w-functions: 


Wo (+9) — (28 + 3 — £27) Wo(g43) + (8 + 1)? w,, = 0 


Soa « SER eee aad (xxix), 
# 1 1 doy, 
@o (541) = (8 + 1) wo, + 4H a (xxx), 
d*w, 1\ dw,, 
eRe (w+ a ae 4 2 (Get Ah) erg 0 sat casagie mae r sd slemeceeaeos (xxxi) 


But the fuller treatment of the w-functions must be deferred. 
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This is also true for s=0, as well as any integer value. It follows accordingly 
that while the total area of any @-function from 0 to o is positive, its & is 
negative for values of s>1. In other words the negative parts of w are on the 
whole furthest from the axis. Again the absolute value of h,, decreases as 

i! 
V2s—1 
to its area to the total mean square deviation of the curve. 

Applying these results to the curve of scatter given by (x), ze. 


1 1 6n—11 50n — 57 
Kaine sae 
bn (2°) (2 dn Oni T99n® + Te00n8 


when s increases, or the higher the w-function the less it contributes relative 


_ 1892 —2125n +2700" 
103,680n' 


Oy — ete.) arte eee (xxxiy), 


we have if A be the whole area and & the radius, 


: ey {t- 3 Le ee a 30 6n—11 all 50n — 57 
Wey  Webyoe Ine! nv 1280 nr? 


2Qr7ro0 2 


77 =—-1892—2125n+270n? ‘ 
= OTRO A ee (xxxv), 


Titegas hee Ue ieee ere 117500 57 
LG ye Walaor MOR ape 3840 n' 


7 1892 —2125n+ 270” ; 
Ug Ge 1 ae en A ete. nee ee (xxxvi). 


: i) : : 
Hence if we even neglect terms of order 7a? We see that the Rayleigh solution 


gives too large an area for the curve of section and too small a swing-radius ; 
these values are 


Rayleigh area, esas Rayleigh swing-radius, o, 
2/270 


i) N etl : : 1 i) 
True area to -, ——/{1— ) : True swing-radius to -, o ( + —) : 
n /Qro 16%” n 


1 
24 8n 

Accordingly for 1 small the graph of the Rayleigh solution tends to exaggerate 
the concentration, 7.e. using it as an approximation we shall somewhat reduce the 
extreme parts of the curve at the expense of exaggerating those near the centre 
of dispersion. 

While there is no difficulty about determining the curve of distribution when 
n is large from (xxxiv), beyond the great labour of dealing with hitherto untabled 
functions, the investigation becomes very troublesome when 7 is small. The 
functions » are suited in this case to represent the discontinuous functions which 
actually form the values of ¢,(7*), but the extreme discontinuity of ¢, (7°) for n 
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small, compels us to use a very great number of -functions, and the convergency 
of (xxxiv) is then small. 

Another method of determining the distribution of the dispersed population has 
then to be applied to the case of small. 


(6) Graphical Solution of the Fundamental Problem for n small. 


Let us consider the general functional relation (i) 


Qar 
Dian 7) = Ey | Pn (re ata (es = 2rl cos 0) dé. 


27 Jo 


Suppose the graph of ¢, from 0 to n/ known. This may be any discontinuous function. 
From nl to «, it will be zero. Let ABD be the graph of $, and OA the axis. 


A 


\B 


ss Oo 


OP=r. Round P describe a circle of radius 7, take the radius PQ, so that the 
angle OP@=86; then clearly, OQ?=7°+l—2rl cos 0; rotate OQ round O down into 
line OD, as ON; draw the ordinate of the graph Ngq, then we have 


Nq = ¢, (7° +P — 271 cos 6) 
and Open Ole = a | "Ng do. 
0 


Hence if we divide the circle up into a number of equal parts, and determine the 
ordinates Vq, corresponding to each of them, we can plot a curve to the base 27, 
of which the mean ordinate will be ¢,,,(O/”), or the ordinate at r of the new 
curve of dispersion for n+1 flights. This can be done for a series of values of 
r from 0 to n+1/ and thus 4,,,(7°) will be determined as a new graph. The area 
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of the plotted curve which gives any new ordinate can be found mechanically. 
It will be seen that the process is theoretically straightforward, but very laborious. 
Thus for the dispersion curve after the fourth flight some 43 points had to be 
found, and this involved the construction of 43 subsidiary curves and their integration. 

There were, of course, graphical difficulties in the construction of the subsidiary 
curves in the neighbourhood of the asymptotes and various devices had to be used, but 
at almost every point there were tests of the accuracy of the work. Some of these 
I shall now notice. 

Case (i). The solution for two flights is: 


Lod en 
$. (7) = rr /4P— 7 ae | Seton Se ee (XxXxXvli). 
=0 r>2l | 


The reader will find no difficulty in deducing this directly from the case of n=1, 
which corresponds to a narrow zone of radius r=/, the rest of the plane being 
unoccupied. Thus: 


b= 5 -from r=l—te tor=l+te ) 
5 


=0 from r=0 tol—te and r=1+4e to 
e being taken indefinitely small. 
By distributing each element of ¢, on the zone round a circle of radius / we 
obtain (xxxvii). 


oe (xxxvui bis), 


The result may be obtained also from (iii) by putting n=2, we. 
Ae ee i ud, (ur) (J, (ul)? du, 
Qa 0 
_ N[(2l+7) x (Ql—-r) ry? 


On 721 (—4) 
=0 from r=2/ to #, 


from r=0 to 21, 


from a theorem of de Sonin by putting a=r, b=c=l. Compare Gray and 
Mathews, p. 239, Ex. 52. 

Case (ii). The solution for three flights may be obtained from that for two, 
by distributing analytically the density given by ¢, round circles of radius / about 
each point. The resulting double integral is then expressible in elliptic integrals™. 
We find: 

: N- -Al T 
RON tes TF KF, K) 
where «= 16l'r/{(r +1) (31—r)}, 
r>Oand <1; 
Nees bt Pi eet Re) eee See (xxxvili). 
Ql Srl @ ) 
where «= (r +1)’ (31—r)/(16?r), 
r>land <3l; 
= 0 r>3ltor=o 


* This solution, or its equivalent, was first sent me by Mr Geoffrey T. Bennett. 
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We have here at r=/ a typical instance of the discontinuity. 


In Table I. columns (i) and (ii) the calculated ordinates of ¢, and ¢, are given, 
the latter having been determined by the use of Legendre’s Tables of the Elliptic 
Integral F. In these cases, as in the later values of the ordinates of the dispersion 
curves, NV is taken as unity. The dispersion curves are plotted in Diagrams I. 
and II. The Rayleigh solution is given in broken line; it will be noticed how 
very far it is from representing the facts at this early stage of the number of 
flights. One of the most interesting features of the investigation is to mark the 
gradual approximation of the discontinuous series of functions to the Gaussian normal 
curve of errors as the value of 7 increases. 


The first test of the graphical method of dealing with the problem was to 
start from the curve for n=2 and construct the graph of ¢, The result was 
found to be extremely close to the elliptic integral solution obtained by analysis 
and calculated from Legendre, and this gave us every confidence in the correctness 
within reasonable limits of the graphical solution, where no such direct verification 
was possible. After the ordinates of any graph had been found their differences 
were plotted, and these difference curves submitted to most careful inspection. 
Larger irregularities led to a reinvestigation of the points, smaller irregularities were 
smoothed with the spline, and from the final smoothed difference curve the ordinates 
were corrected. 


Another test was now possible. In every case 27] 4, (7°) rdr ought to be unity. 
0 


Each ordinate was now multiplied by its 7 and a quadrature formula used to find the 
integral. The integral would usually differ very slightly from unity. Its reciprocal 
was then used as a factor to each ordinate and the ordinates so modified were 
the final corrected ordinates of the corresponding graph. The graphs were made 
on a large scale, and the accompanying Table I, columns (iii)—(vi), gives the 
ordinates of the dispersion curves from four to seven flights. 


Additional tests were as follows: 


Since bral) = | “gba (+P—2lr cos 8) dB, 


Qar 
Bless ke dines (0) = | $, (I) d0 =o, (P), 
0 


or: The axial ordinate of the n+ 1th dispersion curve is the ordinate at a distance J, 
or a flight, from the centre of the nth dispersion curve. Table IV. illustrates 
the degree of accuracy reached here. 


The ordinate at r=/ of the seventh curve given by the expansion in w-functions is 
‘0375, and this is precisely the value of the central ordinate of the eighth curve given 
by the same expansion. Thus the graphical method runs with surprising accuracy 
into the analytical. The Rayleigh solution gives 0398 for the central ordinate of 
the eighth curve as against the ‘0375 of the w-expansion, or the ‘0378 of the 
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TasieE II. Values of the w-functions. 
rio o7 wo o7 Wy o7 wy 77 We a We o7 Wyo C7 Wp 
0 +°159,1550 | +°159,1550 | +°318,3100 | +°954,9300 | + 3-819,7200 | + 19-098,6000 | +114:591,6000 
“il +°158,3611 | +-157,5693 | + -313,5589 | + -935,9497 | + 3-724,9378 | + 18-530,6202 | +110-620,7235 
2 + °156,00385 | +°152,8834 | +:299,5891 | + -880,4201 | +3:449,0302 | + 16:885,5699 |+ 99:177,8011 
3 +°152,1517 | +°145,3049 | + -277,2242 | +°792,4264 | +3:016,3079 | +14:332,2150 |+ 81-601,7164 
“4 + 146,9185 | +°135,1650 | + °247,7634 | +°678,3356 | +2:464,2124 | + 11:127,4046 |+ 59-905,9470 
D +:140,4537 | +:122,8970 | +:212,8751 | +:546,1784 | +1:839,0999 | + 7:583,1582 |+ 36:489,3471 
6 + °132,9374 | +:109,0087 | +:174,4670 | +:404,8965 | +1:191,1906 | + 4:027,9574 |+ 13°802,7339 
of +°124,5713 | +:094,0513 | + -134,5401 | +-263,5330 | + :569,8039 | + -:767,7288 |— 5:975,6743 
8 +°115,5702 | +:078,5877 | +:095,0449 | + °130,4593 | + :016,0640 | — 1:947,9139 |-— 21:205,3207 
oO + 106,1526 | + -063,1608 | + :057,7497 | +-012,7166 | — :435,8815 | — 3:949,8656 |— 30°951,7904 
1:0 | +:096,5323 | + :048,2662 | + :024,1331 | —:084,4658 | — -766,2251 | — 5-161,4614 |- 35:039,7166 
1:2 + °077,4690 | + :021,6913 | —:028,0128 | —:206,6600 | —1:045,7098 | — 5°351,9170 |— 28°874,9613 
14 | +:059,7326 | + -001,1947 °057,3194 | —-235,2026 | — -900,0451 | — 3°455,1198 |-— 12°119,1731 
16 + °044,2510 | — :012,3903 065,5623 | —:194,8306 | — :521,5103 | — -916,7705 |+ 4:126,7483 
1:8 + 031,4966 | —-019,5279 058.4451 | —-119,4328 | — :116,5429 | + 1:050,8391 |+ 12°770,4428 
20 | +:021,5393 021,53938 043,0786 | — -043,0786 | + -172,3144 | + 1°895,4584 |+ 12°751,2656 
2°2 + °014,1523 020,0963 025,8081 | +:013,8001 | + :295,4776 | + 1:723,4411 |+ 7:-400,1858 
2-4 + :008,9341 016,7961 010,9496 | + -043,9712 | + -279,7081 | + 1:008,2741 |+ 1:194,4811 
2°6 + :005,4188 012,8967 | —-000,5180 | +°050,7478 | + ‘188,3692 | + :246,6709 |- 2-829,6050 
2°8 + 003, 1578 009,2208 | + -005,3253 | +-042,6344 | + :083,3863 | — °258,5489 |-— 3-915,1831 
3:0 + :001,7680 006,1880 | + :007,5140 | +:028,5090 | + -003,6465 | — :439,7348 |-— 2°949,4384 
3°2 + °000,9511 003,9185 | +:007,3562 | + :014,7914 | — -038,3979 | — :385,6461 |—  1:308,1481 
34 | + :000,4916 002,38498 | + :006,0410 | +-004,6874 | — -048,6501 | — :231,6523 |+ 007,0283 


I owe this preliminary table of w-functions to the kindness of Dr Alice Lee. 


Much more 


elaborate tables will have to be calculated, if as I anticipate the w-functions are found valuable for 


other purposes. The present table suffices to indicate their general numerical character, and enables 


one to calculate some of the quantities needed in the present memoir. 


graphical construction. The fact that the central ordinate of the sixth curve is 
almost identical with the ordinate of the fourth curve at 7=/, seems conclusive as to 
the general accuracy of the process. 


The above test of the general accuracy of Mr Blakeman’s graphical work is only a 
part of the still more sufficient test that in the seventh curve the graph and the 
w-expansion practically coincide. See Diagram VI. After r=5/ the two curves 
cannot be distinguished, and between r=0 and 3/ the deviation is probably as much 
due to the neglect of higher w-functions as to errors in the graphical treatment. 


Another method adopted by Mr Blakeman for testing the accuracy of his 
graphical work, especially at the end of the range, was to obtain expansions to 
df, (7°), when 7 does not differ much from nl, =nl—€, say, where € is supposed small. 
If f, (€) =¢, ((nl— €)’), then generally for € small : 


a N Eyes ' 
Tn (G)= P (/2)" =a ; ) Oe Eee ora Be BSP ts (xxxix), 


fe | cos! 6d =F (4 (n+ 1)) P(3)/T (3 (n +2). 


where 
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TaBLE III. Table of the v constants. 


m=1 n=6 n=7 n=8 
*) — 041,666,667 | —-035,714,286 | —-031,250,000 
Vg — ‘003,086,420 — ‘002,267,574 — ‘001,736,111 
Ve -000,602,816 000,470,724 000,376,383 
rio 000,104,167 000,067,796 000,046,522 
vi | +000,001,412 | —-000,000,142 | —-000,000,639 
Table of the N constants. 
m=5 n=6 a n=8 
N, — 008,333,333 | —-007,142,857 | —-006,250,000 
N, — 000,123,457 | —-000,090,703 | —-000,069,444 
N, 000,032,600 000,024,174 000,018,636 
Nig 000,000,990 000,000,627 000,000,422 
N,, | —-000,000,078 | —-000,000,047 | —-000,000,033 
m=10 n=6 n=T7 n=8 
N, — -004,166,667 | —-003,571,429 | —-003,125,000 
NV, - | —-000,030,864 | —-000,022676 | —-000,017,361 
Nz 000,008,415 000,006,211 000,004,771 
Ny 000,000,126 000,000,082 000,000,053 
N,, | —-000,000,011 | —-000,000,007 | —-000,000,006 
m= 20 n=6 TT n=8 
NW, | —-002,083,333 | —-001,785,714 | —-001,562,500 
N; ~-000,007,716 | —-000,005,669 | —-000,004,340 
NV, 000,002,137 000,001,574 000,001,207 
Bi 000,000,016 000,000,010 000,000,007 
Ni | —-000,000,0014 | —-000,000,0009 | — -000,000,0006 
TABLE IV. Central ordinates and ordinates at r=l. 


No. of Flights 


Central ordinate Ordinate at r=1 

0 oO 

oe) 70585 
"0585 ioe) 

00 0537 
0537 0537 
0538 0415 
0415 0378 
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a NG as 
~ + (f 2r)"* Bln 2° T (4 (n— 1) 
This can be proved by induction. 


Hence: FS) 


bee tpc pi ten eee RTE (xxxix bis). 


For most of the cases more approximate formulae still were deduced. Thus: 


Fs ie (1 oe os) PME bs ne aa Gin 
10a eal (ee eet (xl) 
HO = ees (F) (14 BF + ee) ores (xii), 
£(9= —Jipl4) (1 ae: i) OP Le 0% eal 


after which the first term only as given by (xxxix) is sufficient. It will be observed 
that after ¢, (7), the curve touches at r=nl or €=0, and the contact becomes higher 
and higher as n increases. Thus, although short of n=, there is no real asymptot- 
ing to the axis, still ¢, (7°) for n> 5 not only vanishes for r=nl, but has increasingly 
higher contact as 7 increases. ‘This explains how the Gaussian curve can fairly well 
represent the state of affairs towards the end of the dispersal range, if n is> 5. 

Mr Blakeman found that the ends of the range for the various cases ran closely 
into the curves (xl) to (xlii), and they were tested and, if needful, corrected by 
these formulae. 

Thus the whole graphical work was kept in check, and, I think, we may be 
confident that the true forms of the dispersal curves for n=4 to 7 are really 
given by our diagrams and tables. 


(7) We may note a few features of these curves. 


Dispersal Curve for Two Flights (Diagram I.). 


There is no discontinuity in the solution from r=9 to 2/, the range within 
which all individuals fall. The curve asymptotes to the vertical at the axis and 
at r=2l. Of course, while the density becomes infinite, the number on any small 
area near 7=0 or r=2/, is finite. Thus the number between the circles of radii 
7, and 7, 18 


2N in UE == ins =) 

(sin 57 — sin 5) 
If r,=0 and r,=«, where e, is small, the number », within the small circle of radius ¢, 
at the centre of dispersion=Ne/(zl). Ifr,=7,+¢, the number lying on the zone of 
Ne, ( ane 


4 — 
breadth e, is —=/ 1-——, , and this if 7, = 2/—e, 1s »,= ay “At the position of 
aa AP oar NL 
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minimum density r,=2I, and the number on the zone 7, to 7, +€, is v,= Ne,/(al/1/2). 
Hence it follows that the numbers on narrow zones «, €, €, in breadth, of equal areas 


ne? = Ale, = 72 2le,, are given by 


Ne|(ml), Nve/(rVJl), and Ne w/2/(xl), 


or in the ratio 


Ne,/(al) :4.Ne,/(al), 4Ne,/(xl) x a 


Thus the total population on a small area at the centre of dispersion is twice that on 
an equal area at the periphery of the distribution, and at both indefinitely greater 
than on an equal belt at the distance of minimum density. The same point can 
be indicated in another way. From r=0 to r=4l is yy of the total area occupied 
after dispersion, it contains ‘16 or about 4 of the total population; from 7=#/ to 
r=2l is 7% of the total area, it contains ‘54.N. In other words the half of the 
area nearest and farthest from the centre of dispersion contains 75 of the dispersed 
population; the “middle” half of the area contains only 38% of the population. 
The nature of the distribution is thus extremely different from that given by the 
rotation of the Gaussian curve about its axis for this small number of flights. 
For in the Gaussian case if the central area we?=277,€,, the area of the zone at 
distance 7,, the population on the centre patch is $Ne/o* and on the zone is 


4Ne?/o° x eT Ho", 
which is always less and diminishes continuously with increase of 7, Thus the 
Rayleigh solution fails in this, as in the next three cases, not only to give the 


form of the curve at dispersion, but to indicate that the dispersed populations on 
zones of equal area round the centre do not decrease uniformly in number. 


Dispersal Curve for Three Flights (Diagram II.). 


The solution is discontinuous at r=/. The density is here infinite, but has 
become finite at the origin. There is no discontinuity at r=2/, but at the end 
of the range the density drops suddenly from a finite value to zero. Thus the in- 
tegral of the Bessel function product (see Eqn. (iii)) is discontinuous at two points. 
The Rayleigh solution is still widely divergent from the true curve of dispersal. 


Dispersal Curve for Four Flights (Diagram III.). 


_ By the rule already referred to (p. 18) the infinite density has returned to 
the origin. There are only two points of discontinuity, ze, at r=land r=4l 
the end of the range, at both of which there is an abrupt change in the slope of 
the curve. The density at the end of the range is now zero and will remain so, 
but the dispersal curve rises at right angles to the axis. The true dispersal curve is 
bending round somewhat to the Rayleigh curve, but the latter is not even yet a 
rough approximation to the facts. 
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Dispersal Curve for Five Flights (Diagram IV.). 


All infinite densities have now finally disappeared. The density vanishes at 
the end of the range, but the dispersal curve makes a finite angle with the horizontal 
axis. There is a marked discontinuity of slope at r=/; a still more noteworthy 
feature is that from r=0 to r=/ the graphical construction, however carefully 
reinvestigated, did not permit of our considering the curve to be anything but 
a straight line. If this could be verified from the analytical expression 


u(r) = | wo (ur) (I, (ult de 


by showing that the integral is independent of r from 0 to J it would be of much 
interest. Even if it be not absolutely true, it exemplifies the extraordinary power 
of such integrals of J products to give extremely close approximations to such simple 
forms as horizontal lines. 

The approach of the Rayleigh curve to the result is now more noticeable. 


Dispersal Curve for Six Flights (Diagram V.). 


There is contact now of the first order at the end of the range. From r=0 
to r=/ the curve of dispersal appears to be a sloping straight line tangential to 
the continuous curve from r=/ to r=6/. No other discontinuity of a low order 
is now visible. The curve, except for the finite slope at r=0, is becoming much 
more of the Gaussian form. It runs fairly closely to the solution in »-functions up 
to w,, in fact is not separable at the extreme part of the range, where the Rayleigh 
curve still gives finite ordinates beyond the possible range. 


Dispersal Curve for Seven Flights (Diagram VL.). 


All sign of discontinuity has gone, the curve is horizontal at the centre of 
dispersion and might be easily mistaken for a normal curve of errors. The expansion 
in w-functions represents the result within the limits almost of constructional error. 
It was not thought necessary to continue the graphical work beyond this stage. 
We may conclude that: 

The deviation of the Rayleigh solution for seven and more flights from the 
true dispersal curve is practically the same as its deviation from the solution in 
a-functions when five terms of that series are retained. 

This I think completes the full solution of the fundamental problem. The dispersal 
curves for the cases of 2 to 7 flights are given in the Table I. of ordinates and the 
Diagrams I. to VI. For higher values the w-function series gives the solution. This 
solution could be applied to calculate the ordinates of the dispersal curve for fewer 
flights than 6 or 7, but several more w-functions would have to be used and the 
arithmetical work—especially while these functions are as yet untabled*—then 
becomes somewhat severe. 


* Table II. provides a preliminary series of values of the wfunctions. 
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(8) Secondary Migration Problems, 


Problem I. On one side of a straight line there is supposed to be a uniform 
distribution of habitats; on the other at starting no habitats. To investigate the 
distribution in the unoccupied area after one migration. Each individual is supposed 
to take n-flights to the new habitat. 


Let YY be the straight line and O a point at distance ¢ from it on the unoccupied 
side of it. Let N be the average density per unit of area on the occupied side. Then 
after an n-flight migration, the contribution from P (co-ordinates 7, y) at O will 
be Nrdxdr ¢, (7’), and integrating this all round a circle of radius r from A to C 
within the occupied area, we have for the quantity F,(c) at O 


cos _ ne 
(= 2 [. )r dydr 
= au cos *¢/r d, (7°) r dr. 


Hence Bia ie) = —2N [ fa (7 oere 2 cies —2N i . Paley) OY. cates. (xliv). 


The evaluation of this integral needs a further consideration of the w-functions. 
By (xiv) 


= oe ne 1 0) pra a 2 | yak 
Chae 5 s\ B) ap \ Be , where B= — 2077/7". 
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Transfer the differentiations from 8 to o° and we have: 


na (o*)$ vr ie rental 


me (0) Hai > thse Soh wis Read ca EE a cee Ie (xlv), 


or, all the w-functions can be found by differentiating the first w-function with regard 
to the standard-deviation squared. Then au (xii) we have 


Wo, = 


ne pee 
ale = fl tn o) agape oP a gy to HV ral aap + fo 
cee (xlvi) 
Thus, if we put o’=t: 
dF : a es 
vale) -2 (1 +S» ie SP oc +(- by Vosl? Te + S x i (?+y’) dy. 
i 2 2 1 ” = (2-+y2)/202 1 —$e?/¢? aa 
But : [, wo, (¢+y’) dy= male (®@+y")/2 dy=s—>¢ 4e?/ Jiro 
nt oslo", 
= ° 
¢c 2/2 clo 
Now i eae ae or [" -2” dv, and this integral vanishes when c=. Hence 
inp 2 a 3 dl 8 dé se — $7 
F,(c)= se (140 Seat te a wae +) |e die 


Since £,(c) clearly vanishes for c infinite, it is not needful to introduce a 
constant. 

It remains accordingly to determine the successive differentials of the integral 
with regard to t. Call the integral 7; then, if n=c/o=c/Jt, 
CG  —c?/2t 


J c 
SS == =T—o,= t. 
oF PON ae A e T= 7Cw,/./ 


By (xlv) we know that d‘w,/dt*=(-—1)%e,,/t. Hence differentiating s—1 times 
we have: 


dy ae Oe; (pe Gig pA leas ) (s—2) Re set) 
dt’ vt (Get 7m de? 34 2! Te eS 
(—1)*" ‘we (s —1) (s—2) 1.3 


e-4 (0, (e-1) + (8-1) e624 + i oT Bee epee 


A MATHEMATICAL THEORY OF RANDOM MIGRATION 27 


at ws Tite s—l1 s—1)(s—2). 
Thus t dé = ( — iy) lay a (ctor ia L! 9 ) OW, (59) Ase her 1 . on OW, (5—3) 
EAE) 18.5 oto n tete-] ee (len 
Substituting in (xlvii) we have, if ¥(n)= ae | j pet hy aie ae eiey 


1 (c) =N E (=) =n oe = {o*w, (dv,+ 20, +42», “E 19 P14) + > 40) 


o 
ap ow, (v, a Ver ZV, a TP Vr) + 525 49) 
2 
+ OW, (v, as $v, “ki $V ats 45.) 


+ 0°o, (v, + 205) + oe Vr.) 
aera WRee valet aye aq Bee ale n 15- (1), 


as far as coefficients of the order v,, and functions of order ,,. 


This is the solution in w-functions. Table IIL, p. 21, gives the values of the v’s for 
certain values of n, and Table II., p. 20, is a preliminary table of the w-functions. 
These will enable us to readily find the values of F,,(c). I have done this for the case 
of n=6 and n=7, which will suffice to illustrate the character of these curves. 
(c/o) can be found at once from Tables of the probability integral. It is drawn 
with a broken line in Diagram VII. and is the Rayleigh solution for this case. I term 
F,,(c) an “infiltration curve” of the first order. 


Substituting the values of the v’s from Table III., we have for n=6: 
F,(c)/N=% (<) + © {-026,712,414 (o°w,) + 053,325,539 (o°w,) 
oO oO 


+ 002,114,303 (o°w,) — 001,029, 898 (o°w,) 
— 000,134,978 (o°a,) — ‘000,001,770 (07a) + ae. ty 
and for n=7: 


F,(c)/N=% (<) ft < {022,850,925 (c°w,) + *045,644,347 (oe,) 


+ 001,578,008 (o%w,) — (000,758,570 (o*a,) 
~ -000,084,525 (c°a.) + 000,000,178 (70,9) + +++}. 


The first term ¥(5) is the ogive curve already drawn corresponding to the 
o 
Rayleigh solution. We see at once that the term ow, will not affect the 
fourth place of decimals. 
4—2 
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TaBLE V. Ordinates of Infiltration Curve over straight Boundary. 


+elo n=6 ni N=@ —elo n=6 n=7 N= 
0 5000 5000. 5000 — = — _ 
1 ‘4614 4612 4602 - ‘l 5386 5388 5398 
2 ‘4231 4228 4207 = 2 5769 5772 5793 
3 3854 3850 3821 - °3 6146 6150 6179 
4 3488 3483 3446 - “4 6512 6517 6554 
5 3135 3128 3085 - °5 6865 6872 6915 
6 “2797, 2790 2743 - 6 7203 ‘7210 7257 
i 2478 2469 2420 - 7 7522 ‘7531 7580 
8 2177 2169 2119 - 8 7823 ‘7831 7881 
9 1898 1889 1841 - 9 8102 8111 8159 
1:0 1640 1632 1587 -10 8360 8368 8413 
1:2 1193 1186 1151 —1:2 8807 8814 8849 
1-4 0834 0830 0808 —14 9166 9170 9192 
1-6 0558 0557 0548 -16 9442 9443 9452 
1:8 0356 0356 0359 -1:8 9644 9644 9641 
2:0 0215 0214 0228 — 2:0 9785 9786 9772 
2:2 0121 0124 0139 — 2:2 9879 9876 9861 
2°4 0064 0066 0082 — 2-4 *9936 9934 9918 
2°6 0030 0033 0047 — 2°6 9970 9967 9953 
2°8 0015 0013 0026 — 2:8 9985 9987 9974 
30 00046 00060 00135 — 3:0 99954 99940 "| -99865 
3:2 00012 00020 00069 — 3:2 ‘99988 ‘99980 99931 
34 00000 00004 00034 -3-4 1-:00000 99996 99966 


m= co is used to denote the Rayleigh solution. 


This table suggests some interesting points. The curves for n=6 and n=7 are 
fairly close together, but differ sensibly from the Rayleigh solution, perhaps 
4 or 5 per cent., where the density is at all material. For many practical 
purposes this might be close enough, and we see that for infiltration as distinct 
from dispersal curves, the Rayleigh solution—owing to integration over an area 
—gives fairly close results. The greatest percentage deviations from the Rayleigh 
solution are to be found in the tail. Now no individual can be found beyond 
the range nl from the boundary, and o=J/4nl; thus the maximum range is J2no, 
or, for n=6 and 7, the maximum range is 3°460. and 3°740 respectively. The 
w-function expansion brings this out well. For n=6 at 3°40 there is not one in 
100,000 individuals, while the Rayleigh solution gives 34. For n=7 there are 
still 4 in the 100,000, because we are a little distance still from the limit of 
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the range. The Rayleigh solution continues to give sensible densities beyond 
the range, although they may be sufficiently small to be neglected in practice. 


For rough purposes a first approximation to the infiltration curves may be found 
from the Rayleigh solution, they will err on the side of safety if we are con- 
sidering the effect of a clearance at a considerable distance from the boundary. 
But with the aid of the tables of the w-functions and the v-coefficients, it is 
not difficult to obtain the actual form of the infiltration curves as I have done 
in the present case. Diagram VII. compares the Rayleigh approximation and the 
infiltration curve for n=7. 


Tt will be seen that an infiltration curve of the first order gives not only the 
density of the population after a first migration into cleared or unoccupied area 
across a straight boundary, but also the diminution of density on the populated 
side of the area, when we put ¢ negative, 7.e. it gives both the ‘depopulation’ 
and ‘repopulation. The reduced density at the boundary is 4N, and if we take 
the point where the infiltration curve cuts the vertical through the boundary as 
origin, we see that it is centrally symmetrical; or the loss of population at a 
given distance from the boundary is exactly equal to the gain at the same 
distance on the opposite side of the boundary. 


If we require an infiltration curve of the second order, we must now multiply 
the ordinates of the curve of the first order by (i) the average fertility of the 
species, say mw, and (ii) the survival rate A. If the environment be the same on 
either side of the boundary, and neither » nor A affected by the density of the 
population, then wA may be treated as a constant and the infiltration curves of 
higher orders can be found with moderate ease for simple cases. We thus have 
the distributions after two, three or more migrations accompanied by reproduction 
and death. On the other hand both » and A may be functions of the density 
of the population, and in this case the ordinates of the infiltration curves of the 
second and higher orders can only be determined when the nature of » and A 
is known. On the whole it is probable that the average fertility depending 
on the mating frequency will be highest where the density is greatest, as mating 
opportunities will then be most frequent, but in such cases the survival rate 
A may be lower, as more enemies are likely to be present and the food supply 
is also likely to be less, where the population is densest. Thus pA as a whole 
may not be very different on the depopulated and repopulated sides of the 
boundary. We shall only consider in this memoir cases in which this product 
is (i) supposed constant throughout, or (ii) constant for each migration season ; but 
supposing uniform environment on both sides of the boundary, it is conceivable 
that »A will be correlated with the population density and this will modify the 
basis of the distribution from which the second and later migrations start. 
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(9) Problem II. To investigate the distribution after m migrations from unai- 
formly densely occugied space across a straight boundary into unoccupied space. 


Let the axis of x be taken perpendicular to the boundary and the axis of 
y be the boundary. Let us consider the density at x=c, on the originally 
unoccupied side of the boundary. Then the density at a distance « from the 
boundary is given by (xlvii), or if we write the operator as @,, we have 


ic.) 


Jip. (2)= NO. = | e?” dx=tm, > (li). 


alo 
Here Q, involves only n and o and not «@. 
Now the distance 7 from the point 2, y to the point c, 0 at which we 
want the density after the next migration is given by: 
Pay +(z—e), 
and pA being the fertility-survival factor, we have for the density at c, 


+o. f+o 
| | pAU,?, (7°) dady. 


ea 

Now Pn (7°) = Q.0, = Q, spp an! . 
To mark that this Q, operates only on this part of the expression, write it 
Q/ and suppose it to operate on o’ written for o. After the operations are com- 


plete we can put o’ again =o. Let 


1 a 8 ae 
1, => ONE 
J Qa ie 
Then if pA be constant (see p. 29): 


AN fa é 2 | p!2 
uy = 0.0! | aie », > oo” beet y3/ da xvdy. 


Completing the integration with regard to y we have: 
wah Q,Q/ [er v, se He aile” da. 
Differentiate with regard to c: 
i ft = HS 0,0/ ike Dee (-e BOOM) de 
Integrate by nin and notice that the part between limits vanishes at both 
of them and we have: 


gh AN = dy, 1 ae 
e re OQ. @e ose Bee Ne 
dy, _ nN 1 6-4/0) 
But ag Ting denarts 


} du, ‘ii pA yuk as we a 
hence : Jae oe OO! seal 6 (Gt i 
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This is integrable and gives: 


Ay _ 1 2 |g Qt q!2 
INTO) Ot wee ee 1 Oe 
de —paA QO; Jin Je to*- 


Integrate with regard to c, and remember that u,=0 if c=; thus: 


=pANQ,Q/ ail ee is a oe BP? +07) aos 


=pANQ, Gna dak ten (lii). 


il @ 
a 
oder 
Comparing this with (li) we see that u, differs from , by (a) the introduction of 
the factors Q,/ and wA and (b) the replacement of o in the lower limit by /o?+o”. 
The process can therefore be repeated as often as we please, and we have for 
Um the value: 


@ 


1 fe 
= (uA)""*NQ,Q0/Q;’... to m terms ef e 2” da, 
J/2a e/a 


where Y=oto?+o"7+... to m terms. 
After the operations indicated by the Q’s are completed, we are to put 
Ce Oe ee 
Now it is clear that a differentiation with regard to any o° is precisely the 


same as one with regard to >. We can therefore write for all the Q’s the 
simple expression 
aa Re Gh 
14% (0 7695-0 get .+(- 1) (0°) a (Sy to 
understanding that d/d(>’) operates only on = and that after the operation is 


completed we can put }=/mo. Thus the complete solution is: 


This is true for c positive or negative, 7.e. whether the density be considered 
at a point on the originally occupied or originally unoccupied side of the boundary. 


Up to terms of order 1/n’ we have for the operator the value 
L+m (ug? — 19 +94 — M9 + P29’) 


-. oe (vq ad Quy. a 2v,v,q°) 


at as pli ) vig’, where g stands for o°d/d (2). 
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Now exactly as on p. 26 we may show that: 


=422 7 0870, NNSA (eae 1h ie? —4$¢2/22 


Wee 


Ua aso 
Er 2m’ Nm@ 3m Yo(e—1) (Mm); 


q : 
Qa p>) 


where: Woig1) (7m) = X2is—1) (n ie se ~ » X2(s—2) (%m) + jee) eae) io ate X2(s—2) (mn) 


sS—1)(s— - 
+! Aerie Ne ic Gia) teee 
nn=e/(Vmo), and y,, is defined on p. 10, Equation (xviii). 
Thus 
1 Sooo Uae | aT 1 1 — me m m 
Un = Ny (uA) 3 (1m) =F De Nn! ke & V4 Ws, (7m) a m3 Ve ws, (2m) 
_ ln tem (m—1) v,} 4 Mitr +10 (m= 1) vy, 
ae We (y s)he pt Ws (1m) 


ea (m—1) v.y,+%m (m—1) (m—2)v? bsp (1m)) ee (liv). 


ms 


We see that this expression converges much more rapidly than that for ¢, (7°), 
if m be at all large. 

The result (liv) might have been reached in a different manner. We might 
have supposed the (wA)""Na individuals to have started from any element a 
on the populated side of the boundary and taken mn flights without multiplying 
to their final resting-place. The effect of this would be that o*=4mml’, and that 
in the values of the v’s we must write mn for n. But doing. this gives us 
precisely the coefficients of the ws in (liv). Thus (liv) is deduced directly 
from (xlix). The proof becomes then much shorter, but it is more artificial ; 
the fact that we may suppose all the unborn individuals to scatter from the 
original centre is not so easily realised, and further it does not in the process 
picture what takes place until the final arrangement after the mth breeding cycle 
is attained. In the method I have adopted we see the exact process of each 
breeding multiplication, its increase of the operating factor by an additional Q,, 
and its increase of the square of the standard deviation by an additional o”. 
Lastly the final form (liv) enables us, without recalculating the »’s for each 
breeding cycle, to see very easily the effect in the case of any n-flight species, of 
taking any number of breeding cycles. 

So long as we keep pA constant of course our result for m breeding cycles 
with n flights will be the same as for a simple scattering for mn flights of a 
larger number of individuals. If wA varies, however, we must adopt the method 
indicated in the above proof, and work out each migration successively. The 
same method must be adopted if a patch be rendered permanently sterile, because 
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in such a case mw is not constant for all parts of the integrated area, and we 
cannot suppose the whole final population to scatter from the original centres. 

If we neglect the %, ,... terms in (liv) we have the value which would 
follow from the Rayleigh solution of the fundamental problem, and this can be 
very readily expressed in geometrical terms. For we mark at once that u, and u, 
are in type identical curves. Take w, and stretch it vertically in the uniform ratio 
of pA”™* to 1, and horizontally in the ratio of vm to 1, and it becomes wp. 
In other words the broken line on Diagram VII. represents the approximate solution 
in this case after m migrations provided we read N(pA)™" for NV on the vertical 
scale and S=/mo for o on the horizontal scale. The Table on p. 33 gives the 
chief results. 

The unit of this table is the length / of “flight.” It will be desirable to illustrate 
its application. Any such application can be of course only a suggestion, and on 
this account the above Table has been calculated to only a few places of decimals. 
But such suggestions may not be without value. They will become more than 
suggestions when our knowledge is greater of the migratory habits of different 
species. At present only rough approximations can be made as to the values 
of m and /, and these admittedly are of small weight. 


Illustration I. In captivity I have noted that H. aspersa will live for five 
years. For two years it does not usually lay eggs, and then it will generally, 
but not invariably, reproduce twice in the year. This is of course subject to 
claustral conditions, and while these seem in some cases unfavourable, in others 
they may be advantageous both in matter of longevity and—owing to the constant 
food supply—in number of broods. This snail, as far as my observation goes, 
appears to return to the same shelter after seeking its food. Leaving such 
“flitters” on one side, I think we might look upon thirty to forty yards as a 
maximum “flight” for such a snail and regard seven or eight such flights between 
its egg layings as on the average an exaggeration. 

We might therefore take /=40 yards, n=8, and an average during life of one 
brood a year as being quite possible approximations in the case of some snails. 

This indicates that the progress across a boundary into unoccupied country 
would be such that 1 per cent. of the density at the boundary and, therefore, 
possibly 4 per cent. of the density in the fully-occupied country, would only be 
reached at 2061 yards from the boundary after 100 migrations. In other words, 
such a species would only progress a mile or two at most in a century. Such 
progress would hardly be noted in any studies hitherto made of distribution; the 
limits of a species a hundred years ago were certainly not closely defined to a mile 
or two, even if they have been recently. Of course there are many other ways in 
which a slow moving species can be transported than by its own “ flights,” and 
further no special stress is laid on the above case, but a study of Table VI. shows 
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that the advance of a slow scattering species* may be comparatively small. The 
inference can accordingly be made that the existing boundaries of the geographical 
distribution of certain forms of animal and plant life which are not marked by 
natural barriers, and which do not correspond to obviously changing environ- 
mental conditions, need not after all be associated with subtle physical differences 
which have escaped the observation of the naturalist. The species may be pro- 
gressing into an unoccupied area, but at a rate hardly observable in the time during 
which accurate distribution observations are available. If this view be correct 
we should expect such boundaries with no apparent environmental change in 
the case of species for which we might reasonably predict a small n and /. 


Illustration IT, I have endeavoured to apply the above theory to the im- 
migration of mosquitoes into a cleared area. We will suppose in the present 
treatment that the area bounded by a straight line (some attempt to allow 
for curvature of the boundary will be considered later) has been cleared but is 
not kept sterile to the species. I shall speak of a district as rendered sterile 
to a species when it is made impossible for it to breed there, and kept sterile 
when the breeding possibilities are persistently destroyed. The distinction is an 
important one, especially in the mosquito case. For in the latter case all mos- 
quitoes are immigrants, and in the former case we have not only immigrants, 
but their produce. 

Major Ronald Ross, who has most kindly provided me with information as 
to mosquito habits, makes the following remarks : 

(a) That the number of mosquitoes produced varies roughly (ceteris parzbus) 
as the extent of surface breeding area. 

(b) That the breeding area can be taken as consisting of numerous isolated 
small pools or vessels of water scattered fairly uniformly over the country. 

(c) That the feeding places (houses, stables, birds, etc.) may be taken as 
scattered pretty uniformly between the breeding pools. 

(d) That abundance or scarcity of food can scarcely influence the question 
much, A single man or bird will yield enough food for many mosquitoes, and 
if they starve it is not because the food is not there, but because they cannot 
reach it. They are therefore not likely to be drawn in general by special abun- 
dance of food in any special direction. Wind tends to make mosquitoes “sit 
tight,” rather than allow themselves to be scattered. 


It would thus appear that on the average an “equi-swampous'’ 
the environment and random “flights” of the mosquito will not be very wide 
of the truth. The difficulty is to form some estimate of m and /. On these 
points again Major Ross came to my help, but naturally the statements he made 


d 


condition of 


were with great reservation. 
* Of course any more quickly moving species that depends on this for food would have the same 


boundary, but in its case the boundary would be environmentally defined. ran 
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(a) From egg to egg (ze. from laying of eggs, hatching, larval and pupal 
stages, to laying of eggs again) takes roughly about a fortnight in hot countries 
with most mosquitoes. In England, gnats may have only one generation or two 
in a summer, but in the tropics they may go on breeding throughout the year. 
In cool countries the egg to egg cycle may be prolonged to a month or two. 
In certain very hot and dry countries, breeding may be checked entirely except 
during the rainy season. I have accordingly taken 20, 10 and 5 breedings to 
the year to represent roughly these conditions. 


¢ 


(8) Major Ross distinguishes between “minor vicissitudes,” which an insect 
makes when it hovers round its victim or mate, and “major vicissitudes” which 
it makes when it passes from feeding place to pool for egg laying. These cor- 
respond to my ‘‘flitters” and “flights.” He considers that they go back to 
water every four or five days, so that a “major vicissitude” occurs every two 
days or so. We might therefore take, excluding flitters, the average number 
of flights to be six or seven. Of course this is the roughest approximation, but still 


not an unreasonable estimate of what probably takes place in the mosquito’s life. 


(y) As to the magnitude of / we have less definite data. Mosquitoes of a 
rare kind have been said to have been found two or three miles from their breeding 
place. Major Ross thinks that Anopheles will exceptionally, when no houses 
are near, probably travel 4 mile for their food, or perhaps further, but he supposes 
the average distance scarcely to exceed 4 mile, and it may, as houses and suitable 
pools often abound not more than 50 yards apart, be not greater, perhaps, than 
100 yards. 

I have accordingly taken 100 yards and 500 yards as likely values for J, and 
considering 1 per cent. of the boundary value of the mosquitoes’ density as a 
limit to their existence and 5 per cent. as objectionable, we have the following 
table : 


Taste VII. Distances from the Boundary of a cleared but not sterile area at 
which 1 per cent. and 5 per cent. of the boundary density of Mosquitoes will 
be found in the course of a Year. 


Supposed number of Breeding Cycles in Year 


5 10 20 


Number of Flights............ 6 7 6 ¥ 6 7 


Density 1 per cent. | = 100 998 1078 1411 1524 1995 2155 
: = 500 4990 5390 7055 7620 9975 10775 


” ” 


Density 5 per cent. | = 100 | 759 820 1074 1160 | 1518 1640 
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The distances are all given in yards. 

Thus we see that the least of these distances for 1 per cent. is greater than 
half a mile, or, if an area be cleared but not rendered sterile, we might expect 
within a year the mosquitoes to reappear within half a mile of the boundary, 
and to reach an objectionable frequency even at this distance for most of the 
cases considered. 

As far then as these rough numbers can be taken to indicate the state of 
affairs, it is needful not only to clear an area but to maintain it sterile. The 
clearance radius may be only 4 mile and is hardly likely to exceed a mile, and 
the above results only mark the progress of immigration in the course of one 
year after the clearance. Further the results would be accentuated if the 
boundary were curved or an approximately circular clearance made. 

It does not appear to me that any substantial difference would be made in the 
main result by reducing n to 8 or 4, although some difference would occur if / 
were reduced to 20 or 30 yards. 


(10) Problem III. Yo determine the distribution after m n-flight migrations 
starting with a centre of population Na, 


The previous two problems indicate the nature of the general solutions to which 
I now proceed. I shall adopt the longer process of proof in this first case as being 
the more suggestive. 


By (xii) and (xlvi), calling the operator as before Q,, we have for the distribution 
at X, Y due to a centre at the origin : 


Ni _~1(¥24 Y2)/o2 
o.(X, Y)=52@, (se b(xr4 Pn) ) ihe le Gas: (Iv). 


Hence the distribution at (h, &) after a second migration of n flights is 


+a to Q en Be — AP + (Y— b/ 
sou (hy B= pa] sb (XY) fee XaT. 


Call the Q, in this Q, and write the o* on which it operates o,’; call the Q, in 
ip, (X, Y), Q;, and the o* on which it operates o,’, we have: 


anos hy+(Y—ky? pee! 


i) os? oy 
hn (h, k) = oy orOn0) ie le "oe dX AY. 
The integrations can be performed and give us 
ie —4 (2+) (o2+os ; 
sha (h, t) =" 01.0, Sa Ae ta) Cina seer (Ivi). 


ee 


This only differs from ,¢, (X, Y) by the introduction of o,+o, for o* and of the 
factor pAQ,.. 


« 
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We can accordingly repeat the process as often as we like and we have: 


A\"™"Na L _2(4242)s2 % 
Pie (h, k) = Va ee Q,, Q,., eee Q;. 32 Cane uta tt eee (lvii), 
where VP=o~tovtost... t+ Ons 


and after the operations have been performed we are to put all the o’s equal 
to o or }’=mo*. But no operator Q, affects any o* in any other operator, and 


d d Raped he 
do? = ds" Thus (a?) doy may be put =(o’) (Sy and this makes all the 


operators identical in form and we may write 
f ee Yes 
QQ Qu [LAM Te —MOY g yt (o) aes 
Ae. ds m 
+ set (—1)'% (0%) Tayt 


a? ae eid 


mee Ca rape (ce gis 


d* 
_1\s Ar 
ici etd eb AGH ay 
In this form of the operator we can now write at once fe and call 


the expression @,”. 


P ad? d* 


Thus Q"=1+N.(2') gem No (BY gmt No) aes) 
Foe (= 1M, (2) 5 oF eee (viii), 
where pects yp Vee 
MW 10 


Ne ea Ae Finis YY 


m* mM 


Mvy+m(m—1)v,7,+4m(m—1) (m—2) ve 


N,= mé 


ol eer (lix). 


These values of the ’s rapidly converge and their values are given in Table III. 
on p. 21 of this paper with those of the »v’s for a few values of m and m. As 
we have seen on p. 82, they are the v’s obtained by using values of nm for n. 


We now have the general solution of distribution from a centre: 


= m-1 Na m 1 —3(+h*)/>? 
mPn (h, k) = (pA) On Q; me ents 


= (wA)" "Na (O,+N,0,+N,0,+ ...+NogQee+...) ...(Lx). 


This is absolutely identical with (xii), except that the constants v are replaced 
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by other constants N of known value, and in every w-function we are to replace 
o by mo* or >, that is to say a uniform stretch in the ratio of J/m to 1 is 
given to any surface z=o,, parallel to the axes of w and y. This is denoted by 
writing ©,, for @,,. 

If we confine our attention to the Rayleigh part of the solution—which will be 
more and more nearly exact as m increases, for the V’s rapidly decrease in value— 
then we have 

Eg Vg EE GAN SoTL A8 Rac aeah tol, Me Ae (1x1), 
and we see that every density gradient curve for the successive migrations is to 
be obtained by a stretch from the first migration density curve. 

In general, however, this result is not absolutely true because the different 
components of the true solution are mixed in different proportions, the N’s being 
functions of m. We see, however, that the stretching rule becomes more and 
more accurate, as we increase either the number of flights or the number of 
migrations. 


(11) Problem IV. To find the form of the general solution for the distribution 
into surrounding space after m migrations of any population imitially spread 
uniformly over any given patch with density N. 

The density at h, k, after m migrations due to a centre Ndady, is by (lvii) 
above 


_, Ndady 1-3 (eh +(y—b)}Y2 
(pA)"? — = O" se oe 
Qa > 


To give the patch let « be integrated from v, to v,, where v, and », will 
usually be functions of y, and then let y be integrated from u, to %. 


We find: 
mln (hy k) = (wasn ake | ‘ | : = e He—- M+ G-ONE doedy ......(1xii). 
T uy, J 0% 
This is the general form of the solution when the population spreads from a 
uniform patch into non-sterile surrounding country. 


If on the other hand we want the distribution after m migrations starting 
with a cleared patch, which is not kept sterile, we have 


FRIES EIS Gere (xiii), 


for the whole district would have had a uniform density of (uA)”""NV had there 
been no clearance. Hence 


F,(h, b)=(wA)"NQ” ( 5 | s | Pa erties FO fea). (xiv). 
mv n ms T 5 a 2 


ib Ug V2 1 Soa ae _ py? ao 
Now F(t, k) =NQ/ (-x| | 2 ee ducdy). 
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Hence the rule: If the solution can be found for a single migration, replace o° 
by mo”, and each v by the proper N, multiply by the factor (uA)”™, and the solution 
for m migrations is deduced. 

It will thus be clear that, if the solution can in any case be found for one 
migration fully, we can at once extend it to the case of any number of migrations, 
with constant fertility-survival factor. 


(12) Problem V. To determine the distribution after a first migration into a 
cleared rectangular area, 


Let the area be the rectangle 2a x 2b, and the origin be taken at its centre 
and axes of « and y parallel respectively to the sides 2a and 2b. Then the density 
at any point h, k, after a single migration F,,(h, k) is given by the principle of 
the last problem by 

Fath, Ry = NEES yn ere eee eee (lxv) 
where F,, (h, k) is the distribution from a uniformly occupied rectangular area 
into surrounding unoccupied space. 


+a (+6 
But F,(h, = N [’ iy $n {(w—h) + (y—B) daedy 
(oxi? ,v-2" 
HOS an es des 
eo; po de x |" 78H gy, 
Qa (Or Nha ait 


Let P,(e) stand for the probability integral 


ae | Pes 
J Dar 0 


il +a a ~“ (a-h)/ a 
Then : = | Gel ae ha win 3% Loe E 
NV 2a70 J -a 2m.) = ee 


i (a -h)lo (athe — 32% 3 
<0 , “ | : )e : 
ap. eo a=") 4P,(2"). 

o o 
Thus: 


F,(h, h=NQ,P, {(*=") P) (“*)} ie (=) +P, (?=*)t. Like 


Now consider the differentiation of P, (2) with regard to o°. 
oR 


da UW d 1 ibe ema lige hab apes 
See PP “)h a8 | a ee a(wi]o?) 
a les 2ado /27 J 0 : ss 20/27 bs 
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Writing o*°=¢ as before, we find 


g E ()} a yi eer ey Ce (Ixvii). 
Hence 
pe ae : ne (s—1) ona) + eo Ue—2) @2(¢—8) : : : 
MTB tna g gto) 
t* £ pe (“)! - ; = (Elem ite (ula Weeds eter. Calon 2M, (Ixviii), 


the expression being the same as that on p. 32. 


Now let us write the following for brevity where n=u/c: 


1 


1 178 
AO esta ro € ae (74%, (9) + 6, (7) + 2 + Vag he s-1) (9) + ---); 


pa aia 
L, (7) fe 2 Nay ne bi (2v,+ BY), (y) auc + S025 Yo (6-2) (7) state )) 


L, (7) =; ee er (3.4 6y,tb, (n) a 0D ihe Daas (7) %) ; 
L,(n) => = noike (40, + 10, (1) +2. BoE) (7) +o.) (lxix), 
and so on. All these functions are directly expressible in q-functions as 
OW ps 27. 
aati fey JES pA Ne ONE) ey sate) ee (Ixx). 
dt* 2 J Qe 


Then we have, if 
n=(a—h)/o, n,=(ath)/o, g=(a—h)/o, «=(ath)/o, 
Fi, (h, k) = N[{Po (1) + Po (12)} {Po (1) + Po (e2)} + {Ls (m) + Ls (ms) {Po (€:) + Po (e)} 
+ {Po (m) + Po(.) {LZ (a) +L, (e)} +1P, (a) +P, (e)} (a (m) +L (2) } 
+... + {P, (e) +P, (e)} {Leas (1) + Leer (ma)p Hee] ccecee ee eeee (Ixxi). 


The Z-functions involve the rapidly converging v-coefficients, and the first few 
terms will suffice to get an idea of the distribution. If we retain only the Rayleigh 
terms we find: 

F,,(h, &) = NTL —{Po (m1) + Po (2) } {Po (4) +Po(e)}] «eee. (1xxii), 
which can be ascertained for given values of a, b, h, k and o from the ordinary 


tables of the probability integral. 
6 
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If we make 6 infinite, then P,(e) and P,(¢)=0 for s>0, and. Z,(e) and 
L(e)=0, Pyle) = Bele) —s anda 
F(A, k)=N {1 — Py (m) —Po (12) — LZ, (9:) — Ly (m2) } vee eee ee (lxxii1). 
This could be deduced directly from (xlix) and it represents the first migration 
distribution into an indefinitely long cleared strip or belt. This is a result of 
some interest as it might approximately apply to the migration into a zone cleared 
by a flood or a fire of certain types of animal or vegetable life. 


(13) Problem VI. To determine the distribution after m migrations into a 
cleared but not sterile rectangular area. 


By the general proposition on p. 39 we have only to write }=mo for oa, 
and the N’s for the v’s in the L’s. Let us put 


m =(a—h)/Z=y//m, y=.) /m, ¢/=e//m, @! =e//m. 
Let 


i 1 1.—3 12 , 
die: (n/’) mar Jon N é ai {NV Wb, (/’) + Now, (m ) Sets Vie Wate—1) (n/’) a ° As 


and so forth, then we have for the full solution: 

mln (h, kh) = wd)" [Po (my) + Po (me) } LP o (e’) +Po(e)} 
ag {L (m,') + Ly! (9./)} {P (e/) + Po (e')} ip {Po (7’) APG (n,')} {L, (<’) + L,! (e’)} 
+{P, (6!) +P (er) } (Le (nm!) + Le (ne) 
+ 0+ {P; (6) +26) a Co) ae) eee ee (lxxiv). 


The terms here will very rapidly converge for any fairly large value of m, 
so that for many purposes we may write the solution: . 


mn (h, &) = (wA)" WN {Po (m') + Po (m’)} {Po (a/) + Po (er')} -200 (Ixxv), 
which can be found at once from the usual tables of the probability integral. 


Illustration I. A rectangular patch 2 miles long and 1 mile broad is cleared of 
mosquitoes, but not retained sterile. What would be the central density at the 
end of the year? Suppose 10 breeding cycles with their scatter migrations, each 
of 6 flights, to take place in the year. Then if we take 200 yards as a possible 
round value for the flight we have: 


m=10, n=6, /=200 yds, of =$nl?=120,000 or o=346°41 yds. 


a= 880 yds, b=1760 yds., », =, =2°540, é, =, = 5081; 
S=/100=1095'44 yds., ni! =n! = *803, e/ =«6/=1°607. 
Hence oH; (0, 0) = (wA)4P, (°803) P, (1°607) N, 
or, using Sheppard’s Tables : 
of, (0, 0) =(pA)'4 x *2890 x -4460N, 


= (pA) x 5156. 
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Thus wf 6(0, 0) =(wA)? (1—*5156) N 
= (mA) 48, 

We see accordingly that if the fertility and the death-rate were the same in 
the clearance and in the populated district outside, the density at the centre 
of the cleared patch would at the end of the year be almost 50 per cent. of 
that in uncleared country. It is thus obvious that clearance can be of small 
use, unless it is followed by permanent preservation of sterility. Even if one 
annual clearance were made it is very unlikely—if the actual values of the 
constants are at all near those assumed—that the mosquitoes would not by the 
9th or 10th breeding cycle within the year before the annual clearance was 
repeated have reached a very substantial density even at the centre of the 
patch. We have thus an additional argument in favour of rendering a district 
not only sterile, but keeping it so. In such a case since vy, and », w,, Ww, are 
negative we shall have a density somewhat less than : 

if, (0, 0) =N {1 —4P, (2°540) P, (5°081)} = N (1—"9889) about. 

Thus : if,(0, 0)="01N approximately. 

It follows that in the centre of such a rectangular patch, there would roughly 
be only about 1 mosquito for every 100 in uncleared country. 

But while this shows that such a sterile patch would be a great improvement 


for a denizen at the centre it is well to enquire what happens in such patches 
some way from the centre. I accordingly add the following illustration. 


Illustration II. A square area of one mile side is cleared and kept permanently 
sterile. What will be the density at the centre and a quarter of a mile from 
the centre on the same assumption as before ? 

Here a=b=880 yds. 

At the centre y,=y,=¢,=6_=2'54 and: 

F,(0, 0)=N[1—4 {P, (2°54)}]=N {1 — (9889)"} = 022N ; 
or, we find one mosquito for every fifty in uncleared country. Taking our 
quarter of a mile directly towards one of the boundaries, we have h = 440, 
k=O, and : 
m=1'27, »,=3'81, g=6=2'54. 
Thus: iF, (440, 0)=N[1—{P, (1:27) + Py (3°81)} {2P, (2°54) }] 
= NV {1 —(-3980 + -4999) (9889)} = 112.N. 


Thus at 1 mile from the centre (or from the edge) of the clearance, the 
density is 11 per cent. of that in uncleared country. It may be doubted whether 
this is a sufficient reduction, and, supposing the above assumptions to be any- 
thing like roughly correct, it may be needful to render more than a square 


mile permanently sterile to protect a patch of one square half-mile. 
6—2 
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On the other hand a cleared but not permanently sterile square mile would 
after a year have a density at the same point—i mile from the centre—of: 


6 (440, 0) = (uA) N [1 —{P, (402) + P, (1°205)} {2P, (°803)}] =(wA) 69.N, 
or of 69 per cent. of that in uncleared country. 


Another point seems of some interest. What is the density at the boundary 
after the first migration ? 


At the middle point of the edge it is 
iF, (880, 0)=N[1—{P, (0) + P, (5°08)} {2P, (2°54)}] 
= N (1 —*5000 x ‘9889) 
= ‘506N. 
This is almost the 4.V of an indefinitely long straight boundary. 
At the corner it is 
if, (880, 880) =N[1—{P, (0) + P, (5°08)? ]="75N nearly, 
or, as we should expect, has risen much beyond the 4 value. 
There is no difficulty in tracing the contour lines of the population density 
in this case. 
If we consider a cycle of 10 breedings in a non-sterile patch we have: 
oF, (880, 0) = (uA)'N [1 — {P, (0) + Py (1-607)} {2P, (-808)}] 
‘T42.N (pA)®, 
and £5 (880, 880) = (pA)N [1 —{P, (0) +P, (1°607)}"] 
= ‘801M (pA). 


Il 


Thus if the patch were not sterile, the effect of the clearance would at the 
boundary after the lapse of a year be marked by a 20 to 25 per cent. reduction. 
The illustrations I have given are of course dependent on the values of the 
constants selected. Such constants have at present been little studied, and 
accordingly small weight can be laid on the actual numerical results. But the 
theory appears to indicate useful lines of inquiry, even if its results will of 
course need to be controlled everywhere by local facts. In a general way there 
can be little doubt that a theory like the present will not only lead to a more 
systematic classification of local facts and to fuller observation of the habits of 
local species, but that this knowledge itself will in its turn test the applicability 
of the theory, or suggest the directions in which it may need modification. 


(14) Problem VII. To determine the distribution after a first nugration into 
a clewred circular area. 


Let the radius of the cleared area be a. Then at distance c from the centre, 
inside or outside the circle of radius a, the distribution Ff, (c) is given by: 
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3) Qar 
F,(c)=N | i . Pace arc COS ONTAOAT mien tn. ene Mteen tae ec (Ixxvi) 
N had ES) Py Aires leone 2 
ipa 0 | if se Bie EO Nge @ = OCS Nee nd Ody 


N 2 3e/o° wo (Qe roe) a\m m 
= = 0, er I i é a” / rs = ( (2) cos” 6dédr. 
a 0 0 u 


oS 
Qar 1/2 
Mee I, os Odd SOE ADS ohh eaten | “cos dO 
0 


> ence eee ! 
_ 4 (2s 1) (2s —8) L7_, 2s! 


2s(2s—2)...2 2 (273))" 


if m be even and =2s. 


Hence: 
F N eee le re) De riios oO y?\8 / e\8 1 s 
n(c)=NQ, <=) i s e as {(5) (5) eam OU eee: (Ixxvil). 
2 S 2s+1 5 o ; 
Now | epi le Te ad (=) =| eB e841] 
a oO oO alo 


aa eval (07 Gag LAr oo Soe Cee (Ixxviil). 


M,,,,(a/c) is thus the 2s+1th moment of the ‘tail’ of a normal or Gaussian 
curve of errors (multiplied by /2z) about its axis. Its values have been tabled 
for s=1, 2; 3 and 4. 


Thus we have: 


ee VEE ed Mom EI / Laie : 
F,, (c)=NQ,e?”! 5 (5) (2° wel hi Seeiaee lye dns (Ixxix). 


But it is easy to see that: 
aorgle tet fy 12 y 1 for AS 
M53 (a/o) =2° sie {1+5 +54 oi) to tare (sa) f- 
Accordingly : 
a Appar iets ea ie 1) +/o8s* Lia xi 1 
Pac NGe 8 (5) {1+ is ( tnntaeg (3) fgay (xxx) 


fom fom 


The successive differentiations of this expression with regard to t=o", involved 
in the operator Q,, which are needful if we wish to give the corrections to 
the Rayleigh solution, are straightforward but extremely laborious. We can 
throw the solution into other forms. 

Write: g=te/o, ~_=40'/o’, 
then we have: 
€, 


€,° 
StS) 


F, ()=NQe SS (Ltet 
0 : 


46 KARL PEARSON 


« 
Here | «e~*dx is the incomplete T-function for an integer value of s. This 


can be found fairly easily from the above series, or may be determined from tables 
of the incomplete T'-function which it is hoped may be shortly published. 


Again : bas wy 2t,/2) =8 a ; 


hence we have: 
fe ()=NQe~ | ol; (Dan) ea) Car esa, REG (Ixxxii), 


a very concise form, which does not, however, simplify the calculations. Integrate 
by parts and we have: 


F, (c)=NQe“+98 © 7, (2iVe)} 
0 2 


io) ls beets 
— —(€,+€2) hy ae he 
ae i d (€€,)° (J, (2Ve.6)}: 
ade . fe ce) gt = Js (24,/z) 
But diet ABS aren eg te (i Jz) ; 
On: ips (c) =NQ, et § {« AU SEO) 
0 (in/e,e,)* 


= NQ,e- +98 (/ =) Jo( Don) eles ieee (Ixxxiii). 
0 2 


This is the solution in Bessel’s functions, and inside the cleared area, where 
e, 1s greater than «, would give fairly good results if tables of the higher 
Bessel’s functions for imaginary values of the argument were available. 

We can also express the solution in terms of w-functions as follows: 


_[° (LEV -ariee rar 

Write ie (=| (5 ) e aan 
1 —$0?/o? 1 oye 
B, (=e 5 5). 

Then F,()=NQS = T, (a) EB, (c) 
0 Ss! 
Now B, (py ae ae =) =270°Sb 5 
s! 2.07 0 ia 


the b’s being undetermined constants, for dividing by the exponential factor we 
have an integer algebraic expression in 7°/o* on both sides. Multiply both sides 
by x»7dr and integrate between 0 and «, p being = or <s. Then: 


1 ) ebay 5 3) rdr = 210°b,, |, Xop pT, 


s! 0 
co rst dr 


=270°b,,(p!)*, by (xix) and (xxi). 
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Therefore by (xvi): 
t —@ 
ow = Sip (P— 1 —s) (p-—2-s) oe (-s) (- ee i B-*el8 dB 


tere {, ave"da, if c= —1/B, 


ee poe). -(s—p+1) 
=(-1) (ply 


Thus E, (rv) = 2:0 { — Sw,+ ayn 1) ea Or. 4 ...(Ixxxiv) 


= 20° U, (r)*, say. 


Now consider: 
= [on BE (= 1) (0 ge, [evra 
=(- (ory gee | ee] 
=(=1)°(0°" Fo: (ao! 


=,, — S@, 


We can now express J,(r) in terms of -functions.. 
We have: 
a4 ee rar 
1,(r)=["s1B, (0) 
s(s-—1) — s(s—1)(s—2) rdr 
=s! 2ro? ie {o.- S@, + (2)? o,— ( (3!) Ws + oes = 


=s12n0"(s-+1) {oy Se eee hee eel es 2) ia | 


iiots 18h BT dl 
= s! 240? (s+1)V,(r). 
Thus F,(c) = NQ,4n%o*S ((s-+1) U, (6) Ve(a))..seceeeeseseee (Ixxxvi), 
0 
8 s(s—1) s(s—1) (s—2) 
where: U(r) =e — pies t (21 O,—- (3!) O,+..., 
s s (s—1) s(s—1)(s—2) 
CS ee SOUL aia Me sta mee tes 


a result which allows of fairly rapid determination from tables of o*e,,. 


There is, perhaps, less difficulty in this form in allowing for the first term 
or two of the operator Q,, for U,(7) and V,(r) can be at once differentiated with 
regard to o°, but even then the final result has considerable complexity. 


* This result involves the expression of any power of 7° in y-functions, 
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The Rayleigh solution value is easily found by putting Q,=1 in any of the 
forms of (Ixxix), (Ixxx), (Ixxxi), (Ixxxiii) or (Ixxxvi). 

A case of peculiar interest arises when c=0, or we take the density at the 
centre of the clearance. In this instance we have: 


F, (0)= Que 30 


Now 
2\2 a? 2\3 d? s 2\s d? 
Q,= 1+», (0°) doy 7) A (Gia ee V5 (a) d(oyt 
and ca a 270° w,, 
2\s ds —ga?/o? eas 2 dw, r d*~* w, 2\s 
therefore (o°) Ta )=270 [ae ape eal (o*) 


= 210" {(—1)*e, +8 (— 1) ev} 
oe (= 1)? {Xa — 8X- (en }- 
Thus 
F,, (0) Ese ges oe {1 — 24x. + (4 — 8%) Xe + (Ye — 4%) Xs + (Ys — 5M) Xe+ +++} 
= 27a?N (a, — 2v,@, + (v, — 8Y,) @, + (V5 — 45) @ + (Vg — 549) s+ ...)...(LXXXVil). 
We are also able to consider the secondary problem : 


What is the distribution into unoccupied space surrounding a uniformly 
occupied circular area due to a first migration ? 


Let the radius of the area be a and let the density at any distance ¢ be 
F,(c) after the first migration. Then clearly, if all space were uniformly filled, 
we should have wniformity after the first migration, or: 

Ff, (ec) + F, (e)=N, 
hence : Ei (Ch = N= Bii (C) facnents eka. teres (Ixxxviii). 


The solution is thus thrown back on the solution obtained for the previous 
problem. In particular at the centre of the populated area we have: 


F.(0)=N=f,(0) oho a eee (Ixxxix). 


We are thus able to calculate the reduced central density due to a migration 
from the area to the surrounding unoccupied district, ze. the effect on population 
of the spread. outwards of a colony. 


(15) Problem VIII. Indirect solution of the General Problem of the 
Random Walk. 


It may not be without interest to put on record the distribution density 
after n flights in the case of a cleared circular area, if it be expressed in 
Kluyver’s manner by the integral of a Bessel’s function product. 
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We have: 


a {2a 
is (=n |. |, b, (C +7" — 2re cos 6) rdédr, 


and 
F,,(c)=N—F, (c) 
if a (2a ) 
=W {1-3 i i | POON ETT PE) Jy (ul)pauirdbar\ 
ea Os de®) 0 


See 
i | i | dy (wr) To (we) {To (ul) du rar | 


| by Neumann’s Theorem (see p. 6) 
E z | i es { i “Sour urd (ur)| aCe du | 


ree | oes | “d(T, (wr) wr} Jy (ue) du | 


0 
by the theorem cited on p. 7, 


=N E ~ | ; Yo ) le furs, (ur) » (uc) de 


N E at | “td Cie ea) du | 


Or, writing v=au, we have: 


F,(=N[1- | “F(T (v 2) {J (v Ale de rae eon 


This expression is concise. The integral expresses the probability that if an 
individual start from the origin and take (n+1) flights, the first of magnitude 
e and the remaining ” of magnitude /, at random, he will find himself within 
a distance a of his starting point. But there does not seem any convenient 
method of evaluating the integral. Comparing with (Ixxxiii) we have the curious 


identity : 


— Lie “(gis An ON . aC : 
[Ft (v | {J (v =) dy ten ri <) J; (2% or) (xe). 


Write c=/, a=r and n—1 for n, then 


< LAN te ie = (2-41) (nm -1) 3 8 5) ( ee ) 
[20 {y. (v=) } dv=1—Q,-1e ° ee J, | 2% Genie 


where Q,,_, is the operator, 


s 


d? di _d 
1+v,(n—1) agp Pe (@— 1) Gat ee H(— 1M (0 — 1) Spee ee 
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or, by (iv), P, (7), the chance that an individual taking n flights from a centre 
should be found within a distance 7 from that centre, is: 


—wli_o oe ere gs!) 7 fo; 2 i 
P (anf ae 8 (i -) we (2 man )} er (xcii). 
: Peet he vei 
Since bn (7) = 5 71 P a ()}, 
we have here the complete: analytical solution in known functions—z.e. the 
Bessel’s functions with imaginary arguments—of my original problem of the 
random walk. But this formal solution provides no better method for shortly 
determining the dispersal. curves than that already indicated in these pages. 


(16) Problem IX. To find the distribution after m migrations each of n 
Jlights, there being originally a circular clearance which is not kept sterile. 


The solution is found by writing mo°* for o°, putting the N’s for the v’s 
in Q, which becomes Q/”, and multiplying by the factor (uA)”* assumed to 
be constant. This can be done to any of the forms (lxxix)—(lxxxiii), or (Ixxxvi). 
If we write: 

é,=e/m and ¢€,=«/m 


we find: 

> Fale) =v GA see Osea S (ay [eersde) WA ee (xciil), 
or: mtn (C) =v (Ay eee ee le (2i/é, x) A) fs ied mae i (xciv). 
Or again: 


J, (207 €,, g) SGM PEE EOE, (xev), 


vias 


~Fu(o) =N (nay enero (,/ — 2) 
aes (0) = (na) 20." 4ntmt oS (U, ae at a) (s+1)) ane eon 


Of these, I have found the first quite as convenient as any other to obtain 
numerical results from. I shall now illustrate the circular patch formulae. 


Illustration I. A circular patch 4 mile radius is cleared of mosquitoes but 
not kept sterile. To find the density at the centre, at 4 mile from the centre, 
and at the margin after ten breeding cycles. 


We shall suppose as before /=200 yards, n=6, and therefore 


2 5 1 @ ° 
o = 120,000 square yards. einai 3227. 


The second term in Q,” will be of the order 34, of the first and I shall 


240 


neglect it. Accordingly the solution may be taken 
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m&n(C) =e &*®) (uA)™-3N (1 +eé,(1+é) + ~ (1 any, =) 


2h 
ah é, ef 


The successive bracketted terms in é, are 
1°3227, 1°3748, 1°3804, 1°3809 and 1°3809, 


which is equal to e** to our number of decimal places. Hence we may put 


oF (€) = (uA) Ner @*%) | L+& (#0582) +5 (¢* — 0061) 


€, Pad oer ce 

tery (ee — 0005) +S e \ 

= (pA)’Ne~ +) (1 — e+ +9 — -05822, — 003022 — -0001€,) 

= (wA)N {1 —e +e-**) (1 — 0582€, — 0030€" — 0001€,)}. 
At centre 

info (0) = (pA) {e-™ (1)} = (pA) -7 24, 
We can test the accuracy of this result by using Equation (lxxvii) which, 
if we put v,=N,, gives: 
wf (0) =(uA)*Ne-® (1 -2.N,x,4+ ...) 
‘6773 


ae ee 9 — "3207 ee 
and X.= 1—€, = (uA) Ne (1+ ihe 2) 


= (pA)'N°730. 
The agreement is accordingly good enough for practical purposes, and we 


may say that within a year the mosquitoes would at the centre of the patch 
have a density 73 per cent. of what they would have in uncleared country. 


I now consider the density at a quarter of a mile from the centre, €,=°0807, 
and using the above formula we find: 
oF (440) = (pA)N (1 — 07" + 07 x +9958) 
Ay ets) 
or, we see that at a quarter mile, midway between centre and boundary of 
the patch, the density is only 2 per cent. more than at the centre. 
Finally, at the boundary itself, ¢,="3227 =e, 
.F, (880) =(wA)N (1 — 07 + e-™ « 9809) 
= (pA) N-79, 
Thus the cleared patch would within the year have filled up with a population 
of mosquitoes varying in density from 73 per cent. at the centre to about 80 per 


cent. at the boundary, or the clearance without permanent sterility would have 


been quite ineffectual with the assumed values of the constants. 
(2 
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Illustration II. Let us assume precisely the same conditions as in the previous 
illustration, except that the area shall be supposed sterile, and we will consider 
what happens at the end of the first migration. 


At the centre we have by Equation (lxxxvii) : 
F (0) = Ne~® {k — 2v, x, + (vy, — 8%) xs + (Y5 — 45) Xp + +} 


But — 2v,= 083,333, X2 (€) = 1—e,= — 2°227,000, 
v,— 8v,= — 082,407, x, (6)=2—4e,+62= — 494,471, 
v,—4v,= — 005,498, x, (€,) = 6 — 18e, + 9c,’ —€, = 8:031,303, 
Vv, — dV, = 000,082, Xs (&) = 24 — 966, 4+ 726, — 1665 +, 
€, = 3'227, = 84°752,347. 
Hence: Lf, (0) = Ne? (1 — 185,583 + 016,024 — 044,156 + 002,850) 


05 [Ve 


This three per cent. of the density in uncleared area might possibly prove 
a trouble and on our assumptions it may be doubted whether the half-mile 
radius is sufficient. If we take the first term only, we find ‘040N, or four 
per cent., not an important practical difference. 


The introduction of even the first modifying term when c is not zero appears 
to lead to such complexity that I content myself with calculating the approximate 
value given by the Rayleigh solution for distances of 4+ and 4 mile from the 
centre of the clearance. In this case ¢,=3'227, ¢,=°807 and 3°227 respectively 
half-way to and at the boundary. I proceed just as before and deduce the 
following approximate value for ,f,(c), 7.e. 

F, (c) =(wA)N {1 —e-* + et) (1 — 20°97 69e, — 7885 le? — 2355; 
— 02286, — 0016e,°—-0001e,")}. 
Hence 

| F, (440) ="179 (uA), corresponding to e,=*807 
and 

iF, (880) ='709 (uA), corresponding to ¢,=3°227. 

Thus the density at 4 of a mile from the centre of the cleared patch would 
be some 18 per cent. of the density in uncleared country. In other words on 
our assumptions a clearance of one mile diameter, if kept sterile, would hardly 
suffice to keep an area of 4 mile diameter free of mosquitoes. 

Compared with a eae boundary, where the density falls to about one 
half that of uncleared country at the boundary, we see that the bending of 
the boundary has a most marked effect in its neighbourhood, the curvature 
raising the boundary density from about 50 to 71 per cent. of the uncleared 
density. In fact the density is almost equal to the 75 per cent. in the boundary 
angle of a square clearance. 
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The differences between a square and a circular patch inscribed in it are 
noteworthy, indicating the marked influence of the area at the angles. Thus 
at the centre we have only 2 per cent. as against 3 per cent., and at + mile 
from the centre 11 per cent. as against 18 per cent. 

As far as the above numerical investigations are to be looked upon as anything 
but illustrations of the nature of the calculations requisite to apply the theory 
of random migration to the mosquito clearance problem, they must be taken: 


(i) As merely an incentive to further study of the manner in which mosquitoes 
scatter from the breeding ponds. It would seem possible, if difficult, to experimentally 
test this by in some way marking a large number of insects, and determining the 
nature and extent of the flight. 


(ii) As indicating that permanent sterility of the protection belt is almost 
certainly needful. The 4 to 3 per cent. of mosquitoes at the centre of the clearance 
amounting to 6 to 18 per cent. at 4 mile distance may or may not be serious, 
but they certainly would very soon be if they were able to breed. 


(iii) As showing that on the rough numbers taken, that a clearance belt 
of probably } mile round a settlement would be the minimum desirable sterile zone. 
But it is quite possible that, when the requisite constants are better known, it 
will be found that smaller belts will suffice. It is possibly rather an exaggerated 
view to suppose a mosquito to make six random flights of 200 yards between 
breeding spot and breeding spot. But certainly many insects I have noted will 
fly with great rapidity in one flight 50, 100 or 200 yards, and these flights are 
quite distinct from “ flitters.” 


(17) Conclusions. The present memoir suffers of course from all the defects which 
must accompany a first attempt to develop a mathematical theory of phenomena 
which have hitherto not been studied with this development in view. The theory 
itself suggests hypotheses and constants which have never yet been considered. 
How far with a broad average of environment in relation to food supply, breeding 
places, shelter, foes, etc. is the spread of a species random? Are any of the 
geographical limits to plant or insect or animal life non-environmental and in 
course of change? If so, statistical studies of the density gradients of such species 
for a few miles either side of the supposed boundary would form most interesting 
work for biometricians. But, apart from this observational work, a good deal of 
experimental inquiry might be usefully attempted with regard to the constants 
of random scatter or flight in the cases of both seeds and insects. 

On the theoretical side there are many problems left untouched. The present 
memoir has only opened up the outskirts of a very big field. It would be of 
value to investigate the number of terms in the expansion in -functions requisite 
to practically reproduce the graphically constructed density distributions for 
migrations of 8, 4 or 5 flights. Our expansion to 6 terms is hardly close enough 
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for practical work until n=6 or 7. Many other shapes of populated or of cleared 
areas would provide problems of some interest, especially when the spread of the 
colony was limited in one or more directions by environmental barriers, such as sea, 
river or mountain range. The problem of sterile areas has been by no means 
exhausted, for in such cases I have only dealt with a result of the first 
migration, but actually there will be a second and later migrations in which 
not only new immigrants will appear but a portion of the first immigrants will 
be emigrants and again able to breed when they reach uncleared country. Our 
solution thus gives only a minimum limit to the percentages if the immigrants 
do not die at the end of the first breeding cycle. Much interest attaches 
also to cases in which the fertility and the death-rate are correlated with the 
density, z.e. 4A is not to be considered a constant. But in these as in other 
problems which suggest themselves, a further preliminary knowledge of some of 
the ecological constants suggested by the present enquiry would be an extremely 
valuable guide to the direction that research should take. 

On the purely mathematical side the problem of the “random walk” may 
now be considered as fairly completely solved. The distribution curves have been 
determined until they pass into an analytical solution expressed by a new type 
of function. The expansion in these functions shows the limits to the accuracy 
of Lord Rayleigh’s solution of a certain allied problem in the theory of sound. 
But the -functions which have arisen in the enquiry have most interesting 
properties, and have led me to a whole series of allied functions of one and 
two variables which I propose to discuss on another occasion. The expansion 
in o-functions will I venture to think be found ultimately to have considerable 
importance for mathematical physics, especially in the evaluation of certain 
definite integrals which arise there. The possibility of practically carrying out such 
expansions depends on the determination of the successive moments (and products) 
of the original function, a process with which every statistician is now fairly 
familiar. But applied to definite mathematical functions it loses the disadvantage 
with which it is burdened in statistical practice—the high relative probable 
error of very high moments—and becomes closely allied to the process of deter- 
mining the integral of the product of any function and a Legendre’s coefficient 
(or solid harmonic). Should the generalised w-functions prove, as I anticipate, 
of some mathematical interest, it will be another illustration of how the need 
of the applied mathematician has thrust him, almost unawares, into the path 
of a novel functional development. 
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First order infiltration curve across a straight boundary. 
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